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Guild of New Zealand, Wellington, New Zealand (proposed by I. L. 
Thomsen); 

Charles Henry Cotter, 9 Atherton Road, Clayhall, Ilford, Essex (proposed 
by E. W. Turner); 

Fahmy Ibrahim Mikhail, M.Sc., 33 Lancaster Court, London, S.E.27 
(proposed by W. H. McCrea); 

*Sydney Nevil Milford, B.A., Ph.D., Yerkes Observatory, Williams Bay, 
Wisconsin, U.S.A. (proposed by S. Chandrasekhar); 

Felix Bernhard Schmeidler, ‘The Observatories, Madingley Road, Cambridge 
(proposed by R. O. Redman); 

Thomas Whitham, 38 Breamish Street, Newcastle-on-Tyne (proposed by 
F. J. Acfield); and 

Ralph E. Williamson, David Dunlap Observatory, Richmond Hill, Ontario, 
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Alan Pennell Lenham, 43 Newcastle Street, Swindon, Wiitshire (proposed 
by D. W. G. Arthur). 

One hundred and twenty-five presents were announced as having been 

received since the last meeting, including :— 

L. M. Milne-Thomson, Jacobian Elliptic Function Tables (presented by the 

author). 
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ON THE STRUCTURE OF COMETS AND THE 
FORMATION OF TAILS 


R. A. Lyttleton 
(Received 1950 April 11) 


Summary 


The self-gravitation of comets is shown to be negligibly small compared 
with the differential force due to the Sun, and hence the separate particles 
of comets move practically in independent orbits, except possibly for comets 
that recede to several hundred astronomical units. A comet cannot be re- 
garded as possessing any permanent shape. The changes in dimensions 
perpendicular to the general orbital plane are shown to be large. Near peri- 
helion every particle of the comet must cross through this plane. Collisions 
between a proportion of the particles will occur and reduce them to finer dust, 
thereby producing particles sufficiently small for the effect of radiation pressure 
to create a tail. Changes of dimensions in the orbital plane during a revolu- 
tion are far smaller. Any slight differences of period must gradually result 
in particles of the comet being distributed right round the orbit. The 
possible contribution of the process to the acceleration of Encke’s comet is 
briefly discussed. 





1. Introduction.—In an earlier paper * the theory of the formation of comets 
was considered and led to the conclusion that comets consist of loose aggregations 
of widely separated dust particles, and that they have masses of order 1018*3 gm. 
and mean densities of order less than 10-!* gm.cm.-*. The linear dimensions of 
comets as observed are subject to considerable uncertainty, for the measured size 
must necessarily refer to the visible portion, whereas it is likely that some of its 
particles extend at lower space density much beyond this limit. Also the observed 
sizes and shapes are themselves subject to curious extensive changes, not merely 
as the orbit is pursued, but sometimes of an irreversible nature not repeated even 
approximately in subsequent revolutions. An extreme instance is the actual 
break-up of a comet into two or more separate smaller comets. Accordingly 
there is difficulty in defining or attaching permanent meaning to the linear 
dimensions or forms of comets, since we have little or no idea just how far the 
particles extend at a density below that necessary for observable effects. Even 
so the observed values will give some general indications of the lower limit of 
possible sizes. ‘The following values are quoted by Chambers f as observed 
examples of large comas: 

TABLE | 
Radius of coma 
Comet 1811 I g*2 x 10!° cm. 
Halley (1835) 4°5 x 10! cm. 
Encke (1828) maximum 2°5 x 10!° cm. 
minimum 2°4 108 cm. 
Holmes 1892 III 1°2 10! cm. 


The object of the present paper is to show how, guided by these general 
indications of cometary masses and dimensions, we are led to simple dynamical 
reasons for the observed changes of shape and the emission from the comet of 


*R. A. Lyttleton, M.N., 108, 465, 1948. 
+ G. F. Chambers, The Story of the Comets, Oxford, p. 222, 1909. 
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material in the form of sufficiently finely divided dust for radiation pressure of 
sunlight to bring about the formation of tails. 

2. Demonstration of the negligible self-gravitation of comets.—If the main mass 
of a comet resides in a central nucleus, the gravitational attraction at points 
outside this nucleus may be calculated as if the comet were a point mass. On 
the other hand, if the central condensation is not appreciable the same method of 
calculating the force due to the comet can be used at points external to it, but with 
less accurate approximation at internal points. If m, denotes the mass of the 
comet, the force at distance r, from it, in gravitational units, will be m,/r?. But 
if R is the distance of the comet from the Sun, the differential force due to the 
Sun at points distant r, apart in the radial direction (for which it is greatest) is 
2 Mr,/R®, where M denotes the mass of the Sun. ‘The comet’s influence will 
therefore be stronger if r,<(m,/2M)P3R. (1) 


As a numerical example, if m,=10'? gm. and R is measured in astronomical 
units (A.U.), this gives r,<4°5=107R cm. Supposing 10” gm. to be an 
estimate of the mass of Encke’s comet, whose greatest distance is 4-I A.U., we have 
r{aphelion) <1-8 x 10° cm., while at the least distance from the Sun of 0-34 A.v. 
we have r,(perihelion)<1-5 x 107 cm. As a second example, if m,=10'* gm. is 
taken as an estimate for Halley’s comet for which a = 17-95, e = 09673, it is found 
that 7,(aphelion)<3-3 x 10® cm. and r,(perihelion)<5-0 x 10’ cm. Since r, 
involves m, only to the one-third power, these values are not much affected by 
moderate uncertainty in the assumed mass. It follows that only particles at 
central distances within these values would be controlled chiefly by the comet’s 
attraction and this only if the main mass lay within the same distances. But the 
values of ‘Table I for the observed maximum and minimum sizes of Encke’s comet, 
for example, are so much greater than these that the attraction of the comet 
itself must be utterly negligible compared with the differential field of the Sun. 
‘This means that particles at such distances, and indeed those much nearer the 
centre of the comet, must pursue practically independent orbits round the'’Sun. 

This result is of such importance to the subsequent argument that it may be 
worth while further emphasizing its validity by calculating the actual ratio of 
the force of the comet on its outer particles to that of the Sun in the above cases. 
This ratio is of order m,R*/2. Mr’. If for particles at the outer edge of Encke’s 
comet when at aphelion r,=2-5 x10” cm., the ratio is about 3 x 10-7, while 
for r,=2°4 10° cm. at perihelion it is about 210-4. Similarly for Halley’s 
comet, using m,=10'* gm. and considering the ratio of forces at distance 
r,=4°5 x 10” cm. in both cases, then at aphelion the ratio is about 4 x 10-4 and 
at perihelion it is about 2 x 107°. 

To illustrate still further the complete validity of regarding the outer particles 
as moving practically independently of the other parts of the comet we can consider 
the time that would be required for such a particle to fall to the centre of the comet 
if otherwise uninfluenced. Thus for Encke’s comet a particle falling from rest 
at initial distance 2-5 x 10! cm. would describe a line ellipse of semi-axis 
1:25 x 10! cm., and with m,=10'? gm. Kepler’s law gives for the time of fall 
about 2000 years. This makes clear that during the orbital period (about 3-3 years) 
the effect of the mass of the comet would be negligible. 

The orbits of the comets show such variety of sizes and shapes that the 
circumstances of each particular one would have to be considered in applying 

19* 
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the above picture. The conclusion will hold at all parts of the orbit for most 
periodic comets, but for a long-period comet moving in a highly elongated orbit 
the self-attraction of the comet at comparable distances from itself will be far 
more effective at aphelion than at perihelion, the ratio in the two cases depending 
on (1 + e)*/(1 —e)*, which may be quite large when e is near I. 

This is perhaps a suitable point to refer to a mistaken impression that seems 
to be held by many writers on comets concerning the so-called ‘tidal influence’’ 
of the Sun. The usual condition for stability of figure against tidal disruption, 
which is necessarily of similar form to the criterion (1) reached above (viz. 
m,/r? > 2 Mr,/R*), is examined and found to fail by a large factor. It is thence 
maintained that comets must be unstable and therefore disrupted at perihelion 
unless there were additional internal forces to prevent it. ‘The argument however 
is not a valid application of tidal theory, for the comet is not a fluid in the hydro- 
dynamic sense; collisions between its particles are very infrequent and certainly 
not elastic. ‘The idea of pressure implied in the condition for the free surface 
to be an equipotential does not apply for a cloud of particles moving independently 
of each other. ‘This criticism does not mean that disruption or similar effects 
cannot occur, indeed it will be shown that they may, but that the question is not 
one to which tidal ideas can be directly applied. 

Accounts of early theories of the changes of form of comets and the formation 
of tails are given in detail by Tisserand.* ‘These theories were based on the 
assumption that a comet consists of a tenuous gaseous atmosphere of negligible 
mass controlled by the field of a heavy central nucleus but also influenced by the 
Sun and in some theories by rotation of the comet as a whole. ‘The forms of 
cometary figures were then investigated by means of equilibrium tidal theory, 
but agreement with observation appears to have been considered by the authors 
to have been far from satisfactory. For instance, as a result of his researches 
on these lines Bessel regarded the existence of an additional repulsive force as 
essential, while Schiaparelli maintained that comets must consist of ‘special 
material’? on which the Sun exerted less attraction than on ordinary matter. 
These studies were made at a time when little or no evidence was available for 
the masses of comets or of their structure. It now seems clear that, in contrast, 
the masses are negligible as far as action on the comet is concerned, and that it 
is not permissible to treat the material as a gas. 

3. Behaviour of particles of the comet near perihelion —The equation of the 
path of a particle situated at the centre of the comet (supposing this to be a definite 
point) may be written in polar form 


a(l —e*) =r(1+ecos§@), (e<1) (2) 


and this will be referred to as the standard orbit, and its plane as the standard plane. 
Assuming the other particles of the comet to have the same period (this is discussed 
later), each orbit will have the same value of a and the same point as focus. 
A possible adjacent orbit for any other particle can be obtained, without change 
in e, by slightly rotating the major axis ASA’ in the standard plane. A possible 
orbit can also be obtained by a slight rotation of the standard orbit about AA’ 
or LL’ (Fig. 1). It is also possible to obtain a near orbit by slightly changing e. 
Lastly, it would be possible for another particle to move in precisely the standard 
orbit but differ slightly in time of perihelion passage. These cases, and any 


* F, F, Tisserand, Mécanique Céleste, IV, Ch. XV and XVI, 1896. 
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combinations of them, are of course representable by slight changes in the 
appropriate elements. 

These adjacent orbits can be regarded as obtained in two ways. First, by 
changes entirely within the standard plane. Second, by a small rotation of any 
one of these first orbits about any line such as NSN’ (Fig. 1) in the standard plane. 
Part of any orbit obtained in this way will evidently lie above the standard plane 
and part below it. Accordingly, a particle describing such an orbit must cross 
through the standard plane twice during each revolution, namely at N and N’, 
Moreover since NN’ necessarily passes through S, it follows that every particle 
of the comet crosses through the standard plane during the time the comet is 
traversing the perihelion side LAL’ between the ends of the latus rectum. 


es 
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For a comet of period P and eccentricity e close to unity, it is easily shown 
from Kepler’s equation that the time spent on the arc LAL’ is approximately 
+(1 —e)°?P, and is therefore small compared with P. For example, for Halley’s 
comet this is about 0:0034P, which is about three months. For some long-period 
comets this part of the orbit may be described in a few days, though the whole 
period may be several centuries. In a time of this order, at most, the whole 
comet must therefore ‘turn itself inside out’’, as it were; particles above the 
plane on the aphelion side of the orbit crossing below, and vice versa. 

There is however an additional feature of nearly parabolic motion tending to 
shorten still further the time during which this occurs for most of the particles 
of the comet. For if, as might reasonably be expected, the nodes N’ on the 
aphelion side, are distributed more or less uniformly round the orbit, then the 
nodes N on the perihelion portion LAL’ will tend to be strongly concentrated 
near A. This will happen simply because in a highly elongated orbit by far the 
greater part of it lies on the aphelion side. In the same way, if the nodes N’ 
were distributed round the orbit more or less at equal time intervals apart, a 
similar result will hold and a strong concentration of the N-nodes occur near A, 
because the comet moves so much more slowly on the aphelion side. ‘The nature 
of the distribution can easily be calculated. 

If n denotes the mean motion, the angular distribution of the points N can be 
found by expressing (1/27)ndt as a distribution function in terms of 6. This 
leads without difficulty to 

n dt 9 do (1—e*)9*(1 + e*) db 
se (1—ecos?6)? 7’ 


2n 





<0<-. (3) 
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From this result the following values have been calculated for Encke’s and 
Halley’s comets of the proportion between — @ and @ of the total integrated value 
between — 7/2 and 7,2. 
: Tasce II 
Encke e=0°846 Halicy e=0°'967 
proportion 6 proportion 
within +6 within +6 
10 O31 10 0°65 
20 0°55 20 0°88 
30 0°70 30 0°95 
40 o-So 40 0°98 
60 °'90 60 0°99 
In both cases the figures show the strong concentration towards perihelion, and 
comparison of the two sets of values shows how rapidly the degree of concentration 
increases as e approaches unity. 

4. The release of material in form suitable for production of tails._We come now 
to a consequence of the theory that seems to throw considerable light on the 
question of the source of the material going to form the tails of comets. The 
difficulty has been to understand why certain comets should more or less suddenly 
begin to emit material for tail formation and how they are able to repeat the process 
almost indefinitely. The present theory appears to go a long way towards solving 
these problems. 

As already mentioned, the actual extents of comets may much exceed the 
observed extents, and the fact that meteor streams of great width surround a 
large part of the orbits of some comets also suggests that comets have greater 
extents than those observed. Meteor showers lasting more than a week are on 
record, and numerous instances of ones lasting for a period of the order of a day. 
This affords an estimate of the widths of streams and thence roughly the overall 
dimensions of comets transverse to their orbits. In a week the Earth describes 
about 2 x 10!" cm. along its orbit and in a day nearly 3 x 10"! cm. Because of the 
high eccentricity of cometary orbits, in most cases the Earth would cross the width 
of the stream more or less directly, though not necessarily through its central 
regions. It is seen that these figures exceed the maximum observed diameters 
of Table | by factors ranging from 2 to about 20. The transverse dimensions 
25H of comets may accordingly exceed the observed sizes by a factor of at least 2. 

Referring now to Fig. 1 (4), let us consider two of the particles of the comet 
initially at points C and C’ on opposite sides of the standard plane when the 
comet is at the end of the minor axis of its orbit and approaching perihelion. 
The figure can also be interpreted as a vector diagram of the velocities of C and 
C’ and of the centre of the comet B when at perihelion. Hence the relative 
velocity of C and C’ will be, with sufficient accuracy, denoting COC’ by 2, 

2i x velocity of comet at perihelion, 
since 27 turns out in all cases to be small and of order 10°? radian. Since the 
velocities of comets at perihelion are of order at least 50 km. sec. it follows that 
the relative speed of two such particles at, or near, perihelion is about § km. sec.~}. 
The following are estimates for a number of comets : 
Tasce III 
. Estimated Angular Perihelion Relative 
Comnet radius 5H size 21 vel. (cm. sec.~) velocity 
Encke ¢x 10°" 6x10~* 7x 10° 4104 
Halley 10"! 3 x10~ 6 « 10% 2x10! 
Holmes *s x 101! 8x 107° 2°5 x 108 2x 108 
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In this table the radius 5H has been estimated by taking, in round figures, twice 
the observed greatest radius quoted in Table I. In some cases higher velocity 
may occur if this underestimates the radius, but for inner particles the speed 
would be less, and also for particles arrested at the standard plane by material 
already at rest there. 

This brings us to the important feature that the attempted crossing of particles: 
through the standard plane in a part of the orbit concentrated near perihelion 
may be prevented by collisions, at speeds of the above order, with other particles 
crossing towards them from the opposite side. A first estimate of the probability 
of a particle undergoing a collision in this way is provided by the brightnesses of 
comets when at great distances from the Sun, at which stage presumably a comet 
shines mainly by reflected light and is in a fairly quiescent state. According to 
Russell, for Halley’s comet only x;4qq of the whole area within its apparent 
boundary is actually occupied by material, the comet being regarded as a two- 
dimensional object projected on the background. If only this proportion of the 
total mass, say 1018 gm., were involved in collisions through the present cause 
it would nevertheless mean about 3 x 10" gm. On the other hand, once such 
collisions start, some of the material will tend to remain near the standard plane 
and increase the likelihood of further collisions, and also some contraction, on 
the whole, of the comet is to be expected near perihelion, as will be seen later. 
This estimate must thus be regarded as a lower limit to the proportion of mass 
concerned in collisions in the neighbourhood of perihelion. 

Now a speed of 2 x 104 cm. sec. ! represents a high velocity judged by ordinary 
standards-——it would require free fall in the Earth’s field from a height of about 
2 km. to produce it-—and there is no doubt that the resulting impulses on particles 
called into play by such collisions would produce stresses far exceeding the 
strengths of ordinary materials. ‘Thus one of the main results of the process 
would be to pulverize the particles of the comet into much smaller particles. 
‘The proportions of different sizes so produced would probably be roughly in 
accordance with a Poisson distribution, though the parameter would be difficult 
if not impossible to settle theoretically, and there is little doubt that a practically 
continuous distribution of sizes would occur. Now it has long been agreed that 
the tails of comets could be produced by the pressure of sunlight on dust particles 
with dimensions of the same order as the wave-length of the light. It has been 
shown by Schwarzschild *, Nicholson +, Proudman {, and others that for a range 
of sizes of the order of the wave-length the outward repulsive force attains to 
about twenty times the gravitational attraction of the Sun. Thus as a comet 
approaches perihelion it will produce, through these simple dynamical causes, 
particles of dust of all sizes smaller than its original constituent particles, and as 
soon as this happens those in a critical range of size will be violently “blown 
away’ by radiation pressure to form the tail. 

In comets that show little or no tail it must be the case that the mechanism no 
longer functions to any appreciable degree. This might result from a general 
reduction in the extent of the comet perpendicular to the standard plane and a 
consequent lowering of the speeds of collision, which would tend to come about 
from the loss of transverse motion in inelastic collisions, and it might also result 
from diminution in the number of particles of sizes exceeding the critical size. 


* K, Schwarzschild, Sitz. d. Math. Phys., Munich, 31, 293, 1901-2. 
+ J. W. Nicholson, M.N., 70, 544, 1910. 
} J. Proudman, M.N., 73, 535, 1913. 
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For if the comet consisted chiefly of particles already smaller than the critical 
size their further break-up would not produce particles within the critical range, 
and no material suitable for tail formation would be produced. 

A further point worth mentioning is the alleged observation of envelopes of 
material emitted with almost explosive violence, often illustrated by drawings 
of Donati’s comet. It seems more likely, on the present theory, that what is 
actually seen is the apparent transmission of an effect through the comet. 
Collisions will not necessarily occur simultaneously at all points of the area of 
the standard plane occupied by the comet, and the envelopes emerging might 
well be the advancing front in the standard plane of the colliding material coming 
from above and below the plane as it becomes observable through being 
transformed to finer dust. 

It is also possible that the collisions would produce considerable heating 
effects. At speed 2 x 104cm. sec.-! a particle possesses kinetic energy 2 < 10° ergs 
per gm. which converted to thermal energy distributed through the whole material 
would raise its temperature by the order of 100 deg. only. But the impacts 
would involve relatively small areas of contact, and energy lost through collisions 
might initially be communicated to only a small part of the mass involved. An 
action similar to that when flints are rubbed together to produce visible flashes 
must accompany the collisions, but with far higher speeds for the particles. In 
such a way much greater heating than 100 deg. might result for a small 
proportion of the material. 

5. Changes in dimensions of comets 

(i) Perpendicular to the orbital plane.—'Vhe line of nodes NSN’ for most 
particles of a comet with e near to unity will lie near the standard major axis ASA’. 
This means that a high proportion of particles tends to pass through the standard 
plane somewhere near A. ‘Those particles that do not do so but have their nodes 
near L or L’ will be so much less in number that they will be too sparsely dis- 
tributed (away from the standard plane) to give any observable coma, and the 
observable extent of the comet will depend on the distribution of particles close 
to the standard plane. Thus a general idea of the changes in size measured 
perpendicular to the plane may be obtained by supposing all the nodes to lie 
precisely at A. In that case the width of the comet 25H, say, at right angles to 
the plane at a general point P(@) will be related to the width 25H, at the end of 
the minor axis by the equation 


6H *e PQ me —e*)!*sin@ 


—— EE ee (. 
5H, BO I+ecos# ') 


For long-period comets with e very near unity, observations as far from perihelion 
as B are not possible, but if the thickness at the end of the latus rectum, 26H), 
say, is taken as unit, the corresponding formula reduces to 


ar = tan $6. (5) 
Since =o corresponds to A, these formulae give theoretically a contraction 
right down to zero at perihelion, but this obviously arises from the approximation 
involved in taking all the nodes at A, so that the results will not be expected to 
hold very near to A. Even so it may be mentioned that a contraction by a factor 
of 100, as between a point near B and perihelion, has been observed for Encke’s 
comet. 
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Systematic measures of the sizes of comets do not figure prominently in 
.available data, but a series of values for Encke’s comet are quoted by Chambers. * 
In the following table are given the theoretical and observed values of 
—log 1005H/5H,. 

Tas_e IV 
(a) Theoretical (b) Observed 
Radius 


ry (A.U.) log 1000H/dHz (10* miles) 
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Ihe general run of agreement between the two sets of values is as satisfactory as 
could be expected in view of the likely sources of uncertainty. It is not stated 
by Chambers which dimension of the comet was observed nor its relation to the 
orbit, though it may probably have been the greatest. (Many drawings of 
Encke’s comet show it as fan-shaped.) ‘There might be some reduction from 
the possible maximum size owing to the greatest extent not lying perpendicular 
to the line of sight. Also the observed size would depend on the mode of 
observation and on the distribution of particles within the comet at each stage. 
The space density of particles probably decreases outwards and there is no reason 
to suppose that the observed limit of the coma closely corresponds to the actual 
limit of the comet beyond which no particles ever pass, if indeed any such limit 
exists. As the comet changes shape, and with this the space distribution of 
particles within it, the observable limit will not necessarily move with the particles. 
In view of these various considerations, comparison of the observed and calculated 
changes of size can only be regarded as of a general nature. 

(ii) Changes in the orbital plane.—-Still assuming that a comet consists of 
particles moving with identical periods, it is plain that, as the orbital velocity 
changes, the spacing of particles must change in proportion. If we consider two 
particles X and Y on the orbit of the centre of the comet at a distance 5D apart, 
the time taken for X to move into the position formerly occupied by Y is 5D, z, 
where v is the velocity at this part of the orbit. If after a finite interval a pair of 
independently moving particles initially at X and Y have reached points X’ and 
Y’ distance 5D’ apart and have orbital velocity v’, then clearly 6D/v =6D'/v". 
If 5D, is the distance between the particles when in the immediate neighbourhood 
of perihelion, it is readily found that at a general point of the orbit 

sp= (1+ 2e cos + e*) . 
I+e 
By means of this relation the following values for 5D/4D,, at various points of the 
orbit have been computed for Encke’s comet and Halley’s comet, and also the 
limiting values for a strictly parabolic orbit. 
* G. F. Chambers, The Story of the Comets, Oxford, p. 223, 1909. 





Pe (6) 
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TABLE V 
Lengthwise expansion of comet : 5D/5D, 

Encke Halley Parabola 

0'083 O'O17 ° 

0°505 0*500 1/2 

0°710 0°707 1/2 

° 1‘ 1'o 1 


This change of dimension is obviously a purely periodic one involving no tendency 
for the leading and following particles to be lost from the comet nor any tendency 
to produce collisions. 

It may be mentioned in regard to these results that few comets are observable 
at positions much beyond the ends of the minor axis, if that, and hence a change 
by a factor of two in its longitudinal extent is the most that is likely to arise 
observationally from this cause. It is also to be remembered that a change in the 
actual overall dimensions of a distribution of particles will not necessarily produce 
a corresponding change in its observable extent. Indeed the two effects might 
even go in opposite ways, for obviously if a distribution of cometary particles 
were dispersed to occupy an infinite volume it would become entirely unobservable. 

6. The assumption of equality of periods of the particles..-Where long-period 
comets are concerned, as the mechanism of their formation shows, these can 
hold together by self-gravitation when near aphelion, and the common periods 
of their particles may thereby be brought about. But in any event only a single 
sweep round the Sun of a long-period comet is observable, and there seems to 
be no particular feature likely to develop, during the short portion of its arc 
for which observation is possible, from any slight departure from the standard 
period of a proportion of its particles. 

For short-period comets, however, self-gravitation must be always negligible. 
It has long been agreed that these comets have come about through deflections 
by Jupiter of long-period comets. In any such deflection, because of the finite 
size of a comet, a distribution of slightly different periods would be produced 
among its particles. Now the effect of a slightly different period for any particle 
will obviously be to cause it to move gradually out of the comet and either go on 
ahead if the period is less or fall back relative to the comet as a whole if the period 
is greater. In the same way collisions within a comet of its particles may cause 
slight changes in period and again lead to a gradual dispersal along the orbit of 
such particles. It is of course well known that meteor streams occupying the 
whole or a large part of the orbits are associated with certain comets. On this 
picture a comet is therefore to be regarded as a distribution of particles extending, 
in the orbital direction, far beyond its observable limits. ‘The observable portion 
probably corresponds to the small densest region (per unit area) of this distribution. 
Short-period comets are faint objects and probably represent only small pro- 
portions of the original long-period comets from which they have survived, the 
remaining part having gradually become dispersed along the orbit. 

7. The acceleration of Encke’s comet.—The collision of particles coming from 
opposite sides of the standard plane must result in slight changes in the orbital 
elements for the motion of the statistical centre of the mass as a whole. If we 
consider two equal symmetrically placed particles, such as C and C’ of Fig. 1 (4), 
their angular momentum vectors before collision are inclined at angles i and —i 
to the normal to the standard plane, but after an ideal inelastic collision both. 
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would be perpendicular to it. Supposing for simplicity that collisions occur at 
perihelion, the changes will be such that a(1—e) remains unaltered but the 
velocity there and angular momentum are reduced by a factor cost. Actually 
for the comet only a proportion of its mass would be involved in such collisions, 
but the general effect of the process on the orbital elements must nevertheless be 
that of a decrease in h without change in a(1—e). From the usual formulae for 


elliptic motion it is readily found in the present case that 
bn 38a (1 +e) Oh 


n 2a I—-e hk 


be _ bh dl oh 
ee ee 


For Encke’s comet e=0-8458, and the changes are in the following ratios: 
dll dee _daja__dnjn___ oh (8) 
2 4:36 239006CU — 35°99 
‘These values are substantially different from the corresponding changes due to 
the Poynting-Robertson effect *, for which the values are 
dll gee daa bnin 
7. te Hae (9) 
when adjusted to agree for 5///. Robertson himself concludes from his study 
of the effect of radiation pressure that the drag produced leads to values that 
cannot be reconciled simultaneously between the several elements. But the 
present process provides a further mechanism that alters the elements in a different 
way. It seems probable that a suitable combination of the two effects could lead 
to a resolution or improvement of the difficulty referred to by Robertson, but its 
discussion would go beyond the object of the present paper. 
St. John’s College, 
Cambridge : 


1950 April to. 


* H. P. Robertson, M.N., 97, 436, 1937. 
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Summary 


An analysis of rotational currents is made for the case of small oblateness 
and constant angular velocity,w=const. Contrary to von Zeipel’s paradox, the 
average energy generation over a level surface of the rotating star is subject to 
no restrictions; the balancing convective circulation is derived by developing 
Eddington’s method in more detail; a complete system of circulation is 
calculated for a particular model. As compared with the speed of nuclear 
reactions, the mixing efficiency of the rotational currents for w=10~* (or 
30 times the solar rotation) is almost negligible, being of the same order as gas 
diffusion and being capable of mixing efficiently not more than the innermost 
20 per cent of the stellar mass; for the solar angular velocity the effect is 1000 
times smaller, thus practically nil. Eddington’s suggestion that rotational 
currents are capable of maintaining uniform chemical composition in the 
radiative equilibrium layers of a star must be definitely rejected except for 
configurations on the verge of disruption by rotation. 





1. From a paradoxical theorem of von Zeipel (1), Eddington (2) and Vogt (3) 
concluded that in the radiative equilibrium layers of a rotating star, convection 
currents must necessarily arise; these are forced large-scale currents working 
their way against the stable temperature gradient. Since then the question has 
been treated by different authors on several occasions (4, 5, 6, 12). Eddington (7) 
gave a formula to estimate the average velocity of the currents; however, for the 
Sun’s angular velocity of rotation the velocity is far too slow to produce appreciable 
mixing of the material (8). 

In view of the fundamental importance of this question from the standpoint 
of stellar evolution, let us reconsider the von Zeipel problem of constant angular 
velocity of rotation for the case of w small (small oblateness). We start from 
Eddington’s equation (2) (7): 

f@)(47Gp — 20%) —f ()(d/dn)? =ple,+e.), (1) 
where ® is the potential (from combined gravitation and rotation), g = — d@/dn 
is the surface gravity and 

O O 
f(®) = - a (2) 


oO 
5 


is a pure function of ® (g) with QO denoting the outward flux of radiation in the 
direction of the normal per unit of time and area of the level surface; p is the 
density, «,+¢«, is the energy generated per unit of mass and time. Evidently 
f'(®) = df/d® is also a pure function of ®, thus remaining constant over the level 
surface. 

Von Zeipel’s paradox follows, when we assume that «,+ «, is constant over 
the level surface; as (d®/dn)* is not constant, its coefficient in (1) f’(®) must be 
zero, f(®)=const. and from (1) 


ad 
€yt+eg~ (x- &). (3) 


* Received in original form 1949 December 1, 
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This conclusion has been criticized ordinarily from the standpoint of physical 
soundness; the intrinsic energy sources will never follow a law like that, and 
convectional transfer of heat has been advocated to obtain the peculiar prescribed 
form (3) of the energy generation. In this case, of course, convection is treated 
as a supplementary source of energy. However, the only conclusion so far is 
that «, +, cannot be constant over a level surface of a rotating star. 

Let «, be constant over a level surface, ¢, =«€,(®) being a pure function of the 
potential, and «, be the variable component of energy generation, its zero point 
being chosen in such a manner that the mean value of e, over the level surface 
is zero. Thus, over a level surface ¢,;+¢,=¢,. Consider the case of small 
oblateness, when 

wry 


a= ry (4) 


is small (7, and gy being certain average values of radius and gravity for the level 
surface), so that squares and higher powers of « can be neglected. In such a 
case g = —d®/dn involves a second zonal spherical harmonic of the astrocentric 
co-latitude 7/2 —4, 


dD - 9 
~ 5 =8 =8ol1 + Aa(sin® $ — })], 


where 4 is a factor between 2 and 1:5, gg the acceleration of gravity at r=79, and 
sin® dy = 4, 69 =35° 15’; the average over the level surface is evidently 


z= ( J" gsing de) / ( |" singag) ane 


According to (1) and (2), with (5) as the only zonal variable, both eg and Q 
involve second zonal spherical harmonics similar to (5) but with different amplitude 
factors; their mean values over a level surface coincide therefore with their 
true values at r=7p, 6 = dp. 

When = is small, the angle (m,r) between the normal and the radius is of the 
order of x, thus cos (n, r) differs from 1 by a small quantity of second order. ‘Thus 
dd om ee ea 
= cos (n, 7) = > 

may be assumed. 

In (1), (d®/dn)* =g>5 +(g* —gj) consists of a constant part gj and a variable 
g’—gj, the average of which is zero, according to (5) when «* is neglected. 
Separating the terms which are constant from those which are variable over a 
level surface in (I), we obtain two conditions : 


S()(42Gp -. 2w*) —f'() . g5 =pe, (6) 
and 


—f'(®)(g° — 8) = —f'(®) - 8 (f ~ :) = Pep. (7) 
Along the radius at ¢=¢p, the gravitational force is equal to the attraction of a 
spherical body with r =r, and the same mass and radial density distribution; thus 
the potential is 

ae = 
7 dry + 4w*r, cos” do, (8) 


Tr “0 


=G| 
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M denoting the mass inside the level surface of ro, and 
ro =1 (1 —ksin? $o) =r,(1 — 3h), (9) 
where r, is the equatorial radius, k the oblateness (a —), a. 
From (8), with cos* dy = §, 
d® GM . 
~, dr, =25= a a $w"*Fy. (10) 


With Q,=L/4nr;, L denoting the total radiation flux through the level surface, 
(10) inserted into (2) yields 


f(@)= — ‘ 


477g, 4n(GM— Gur)’ 





(11) 


Further, 
(0) = LO) /@ __ + HO) 
0 0 So 4 
or from (11) and (10), with dL. = 47rzpe, dry and dM = 4xr7pdro, 
- 2 Gel wl oe 
J (®)g5 = — pe, t+ 72 te 2nr2gy (12) 

Substituting (12) and (11) into (6) we find that everything cancels, and that (6) is 
satisfied by arbitrary «,. No condition upon ¢,, the principal component of 
energy generation, is imposed by equation (1). The average energy generation 
is not affected by rotation. 

Substituting now (12) into (7), a formula for the convective energy generation 
is obtained, with the aid of (10): 


— €s 


("- src) 
€-1) € €| 2nGp 


w" , 
es sth) 
here «, = L/ M and p= M (4zr;, 3) are the average interior energy generation and 
density. Formula (13) is essentially Eddington’s (7) formula (8) but with a little 
more detail. In the equatorial plane g, <go, and as «, <€, according to our present 
ideas of stellar energy generation, ¢, is positive; energy must be released along 
the equatorial plane which, in a stable stratification such as a radiative equilibrium 
layer normally is, requires a downward motion along the equatorial radii and an 
upward polar current, with a dead line at 35° latitude. 

No condition regarding molecular weight has been imposed; therefore the 
conclusions are valid also in the case of a radial stratification by molecular 
weight. ‘The rotational currents indicated by (13) will work toward equalization 
of the chemical composition even when there is stable stratification by molecular 
weight. However, in an exhausted central portion of an advanced stellar model, 
the hydrogen abundance may increase outwards so fast that €, may increase, in 
spite of the decreasing temperature and pressure ; in such a case €, >-€, and by (13) 
the direction of the current will be reversed—up at the equator, down from the 
poles. 

Following Eddington we may evaluate the vertical velocity of the current; 
as there is no actual excess of temperature, and as the convection is required to 
balance the prescribed amount transported by radiation, no allowance for a 
lateral loss of heat is to be made, as should have been done in the case of thermal 
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convection (13). Let unit mass be moving so as to change its state according 
to a polytropic index n; n is prescribed by the conditions of radiative equilibrium. 
For an infinitesimal change of state dP, dV (V =1/p), dT, the external work is 
equal to the sum of dq, the energy released (transported away by radiation), and 
the kinetic energy of the molecules stored inside, 
—~ PdV =dq+c,dT. (14) 
‘Together with the gas equation PV =(y—1)c,7, where y=c, ¢,, and with the 
polytropic equation Px T”*!, this yields 
dq =c,dT(|n(y—1)—1] -e,dt(* -1), (15) 
where n, =1/(y—1) is the “‘adiabatic”’ value of the polytropic index. Following 
Eddington (7), call 57 the potential change of temperature under adiabatic 
conditions. According to (13), equation (30), 
I 
n, 


whence d7T'= - er = dT); 
a) 
substituting this into (15), and with 1+ 1.”,=y, yc, =c,, we obtain 
dq=c,(5T—dT). 
Eddington (7) uses here c,, for c,, which needs to be corrected. 
Let the vertical velocity of the current be v,=dr/dt. Evidently 
€, =dq/dt =c,(dT/dr).v,|n(y —1)—1], 

by (15). With the aid of (13), equation (24), we obtain 

dr €> €.(n +1) 


=P 


dt" &{n(y—1)-1]  —g(n—n,) 





; (17) 


where £=d7T dr. Eddington (7) obtains v, too large by a factor of y, as the 
consequence of the above-mentioned inaccuracy. 

The effective acceleration of gravity varies over a level surface according to (5), 
with gy given by (10). For a complete concentration of mass in the centre the 
theorem of Clairaut yields, for an oblateness k= 4a, A =2; for the Earth, the 
observed values are k=0-g80x% and A=1:539. A=3—hk/ is a relatively slow 
function of the density concentration; with Chandrasekhar’s (10) calculations 
for polytropes the run of the oblateness ratio k/« is as follows: 

Polytrope indexn ... 4 3 2 15 ro ~—- (Earth) 
Pc!Pm 00 622 54°2 II-40 5°99 3°29 (1°8) 
k/« Or5 0°502 O°515 0°575 0°639 0°760 a-g8o 
Interpolation O-5 =Or501 0-516 0575 0643 0-760 0975 


The last line gives an empirical interpolation according to k/« =0-5 + 0-850p,,/p,; 
the agreement is excellent, even in the case of the non-polytropic Earth. We may 
therefore assume, for a level ellipsoidal surface with p as the mean density inside 
the surface, 


0-856p 
ms 


A=2- 


(18) 
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Although (18) is not meant to apply to an interior level surface, we may safely 
use it also in that case, the approximation being better than, e. g., the assumption 
A = const. 

For a nearly spherical star, « being small and g=GM,/r* practically, the 
expression for the vertical component velocity of the rotational heat exchange 
circulation becomes finally (from (17), (13), (5) and (18)) 


I-712p 


See)? 40 (hls (@-1)]—choies-ns 


the factors containing %« being retained, in order that the formula be valid also 
for the lowest values of the density. 

For an interior point of the Sun when w =3 x 10-*=const., M,=0°3Mo, 
r=02Ro, g=2 10°, «=6°7 erg/gm.sec., €,=0, N=2°5, N,=1'5, 1 =6°3 x 107%, 
pip.=4, we get v,=2x 10°! (sin?d—4) cm./sec. This current velocity is 
negligible; the maximum value in the direction of the pole, 6=90°, would 
correspond to a displacement of only 100 kilometres in 3000 million years. Even 
with a tenfold value of w the mixing effect would still be negligible, giving a 
displacement of 4, of the solar radius in 3 x 10° years. It should be safe to 
conclude, in the case of the Sun at least, that the mixing efficiency of the 
rotational currents in the interior is practically nil. 

For an external subphotospheric region ((14), B)), #=3 x 10°*, r=o-95R 
M,=Mo, g=3x104, €=2, €,=0, m=2:0, p,/p=63, «=2x 10-5, we obtain 
v, =3'2 x 10°8(sin?d—4)*, thus a value 160 times the value for the interior. 
In 3 x 10° years the maximum vertical displacement is 2 x 10° cm., a value which 
is small but not inconsiderable, amounting to two-thirds of the depth of the 
layer below the photosphere; thus, in the peripheral zone of the Sun rotational 
currents may play a noticeable role. 

In fast-rotating stars conditions may be different. Our formulae may be 
used for order-of-magnitude estimates even with large values of «. ‘Take for the 
same interior point as before w=1-5 x 10-* (period of revolution 1-2 hours), 
% =0°25, k~o-25 (oblateness), thus a marginal case near instability; this gives 
v,=6 x 10°5(sin® ¢ — 4), covering one-fifth of the Sun’s radius within ten to twenty 
million years. ‘The mixing efficiency would appear quite considerable but for 
the inevitable stratification which would limit the circulation to separate zones, 
without much exchange of matter between the zones (cf. (g), p. 286, and following 
pages of the present paper). 

2. The Coriolis deflection of thermal convection in a rotating medium breaks. 
up the circulation into small-scale vortices (cellular or irregular) (13). It is 
otherwise with the rotational currents in a thermally stable medium, considered 
in the preceding section. The vertical motion is the same over a whole zone of 
latitude; the medium as a whole is lifted or lowered, supported by the radial 
force of pressure versus gravitation; it is a case of motion under a purely central 
force, in which case the law of conservation of angular momentum 





ws* = const. (20) 


applies, where s=rcos¢ is the distance from the axis of rotation. This applies 
equally to the radial and to the meridional displacements which are the necessary 


* Valid only when p/p< 10'~a-!. If this is not the case, the sign of v, will be reversed. 
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consequences of the radial movements. No vortices will be created except by 
friction which, however, in the absence of turbulence, is small at stellar dimensions. 
Differences in the angular velocity of rotation will be created which change the 
basic picture of von Zeipel’s problem and ultimately influence the result; at 
the start, however, the conclusions of Section 1 should be valid and, as the Sun 
cannot have proceeded far from the start with respect to rotational currents, in 
that particular case we are entitled to apply the theory with some confidence. 
Considering that peripheral turbulence makes the photospheric rotation follow 
the rotation of the underlying radiative equilibrium layer (cf. also (12)), and 
assuming for the latter a radial outward displacement in the polar region of 
2x10 cm. or 4 part of the radius as estimated above, we should expect in 
310° years a polar slowing down of rotation by 4 and an equatorial 


acceleration by + of the average, at photospheric level; this would give a 
difference in the period of rotation of 2:2 days, about 25 per cent of the observed 
difference. ‘The agreement can be improved by setting the top of the radiative 
zone somewhat nearer to the surface than assumed in the present sample 
computation. 

Here only radial displacements are considered; this is correct for the pole 
and the equator but not for intermediate latitudes. However, for an incomplete 
circuit as in the present case, meridional displacements directed from the pole 
toward the equator will always lead to an equatorial acceleration and will not 
change the picture. ‘This apparent paradox may be easily understood if we 
consider that, although the meridional displacement away from the pole leads 
to a decrease in the angular velocity of rotation, the decrease is insignificant at 
the equator and is largest near the pole; thus the resulting angular velocity 
will be largest at the equator. With s=rcos¢é and r=const. (20) yields 
dw=2wtan¢ddd. If the original angular velocity, we, is assumed to be 
uniform, the resulting angular velocity after a frictionless displacement Ad 
iS W=Wy+2watandAd. For Ad<o, or a displacement toward the equator, the 
angular velocity will steadily decrease from equator to pole, for any reasonable 
law of variation of Ad with latitude. 

3. The structure and stratification of the rotational circulation at an initial 
stage can be easily derived from the condition of continuity of flow, for the case 
of small oblateness as before. ‘The radial velocity of the forced flow as given 
by (19) may be written as 

v, =vp(sin?d— }), (21) 
where vp» contains all except the zonal harmonic factor of (19) and is constant 
over a level surface, depending solely upon the average radius r of that surface. 
The radial flow of mass through a zone between ¢, =0 and ¢, =¢ is 

o 
F,= one 2nrcosd.rdh= — inpr*vy sin ¢ cos* ¢. 
d o, 


The total inward flow through the equatorial zone between ¢,=0 and 
$2 =$9 =35 15’ is 


47 
Fy= — —=pr'v, 


9v 3 
being equal and opposite in sign to the total outward flow toward the polar cap 
(between ¢,=¢,) and ¢,=90°). The flow is of the same sign as v,, i.e. to be 
reckoned positive outwards. 


20 
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Let v, denote the tangential meridional component of velocity, positive when 
¢ increases (motion from equator to pole). The condition of continuity requires 
the excess of the radial flow to be removed by the meridional displacement, 
dF, = —2nrcos¢.dr.pv,, or with (22) at ¢=const. (conical lateral surface) 

_ , sngdcos¢ d_ , 
v= + SESE 5, (ore) (24) 
The meridional or lateral velocity through the conical surface ¢ = dy (sin® do = 4) is 
ad v2d 
. mm gpr a (Pr *v9); (25) 
at this surface v, =o, thus the flow is purely lateral. 

The third rotational component at right angles to v, and v, does not influence 
these results as long as formulae (21) or (19) remain valid. 

With (25), (24) may be written in the form 


v,= Jom sin d cos d, (26) 
V2 


v,, being a constant for a given level surface. 
The equation of a stream-line (projected rotationally on a meridian plane) 
follows then from dr/dt =v,, rdé/dt =v,, (21) and (26) as 
A(r) — (ro) = A(¢), (27) 
where 


6(r) = (24 al 


m (28) 
with 

A = —0:3618[0-4147 + log (sin ¢ cos* 4)], 
r=1o, $=» being the starting coordinates for the integration; ry is a maximum 
or a minimum radius of the stream-line. 

For a certain composite model “ B,"’ of uniform composition (14) the stream- 
lines and the time intervals (t)) covered by the circuit are calculated. Table I 
contains the basic functions (the complete description of the model is given in 
another paper (14)), Table II the characteristics of the stream-lines, all 
calculated with w= 10~*, corresponding to a period of rotation of 17-4 hours or 
an equatorial velocity of rotation of 63 km./sec. and an oblateness of 0-00912 (at 
photospheric level). The boundary condition chosen is logp,=8-572, at 
T.,, =6406°, being close enough to the photospheric conditions of a model solar 
atmosphere i in which absorption by the negative ion of hydrogen is taken into 
account (11).* 

The resulting system of rotational circulation is represented graphically in 
Fig. I. 

There are two major circulations; one, linked with the core, limited by 
r = 49°32 x 10° cm., which contains 99 per cent of the mass, its stream-lines 
emanating from the pole of the core and returning to its equator; an isolated 
circulation of not much importance is included, around a dead point at ¢ = 35°-2, 
r=29°7 x10. Another circulation, connected with the peripheral convective 
regions, occupies the zone outside r = 49°32 x 10°. The stream-lines meet the 
equator at different angles from 0° to go°, whereas they leave the pole tangentially. 
At the boundaries of the core (r =7-58 x 10°) and of the peripheral convective 


* For 3 per cent metals Miinch’s data yield log Pp =8-70 at T,=5713°. 
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region (r=6-1 x 10") singularities arise, as the velocity components become 
infinite (n=n,, cf. equations (19) and (24)); this means that our equations of 
static equilibrium become invalid near these critical surfaces, as soon as v, or 
v, cease to be small compared with wr; this happens, however, in a layer some 
metres thick and is of no practical importance. Enhanced singularities appear 
at the pole and the equator of the core (and of the bottom of the peripheral 
convective zone) which, however, are equally of little consequence. 








Fic. 1.—System of rotational circulation in Model B),. 
The figure represents a meridional quadrant with the stream-lines projected thereon. 
©: Dead points (dead rings). 
C: Central convective region. 
S: Sub-photospheric convective region. 
E: Equator. 
P: Pole. 

More important than the spatial stratification is the “ stratification in time” ; 
the time of circulation varies from 107 to 10! years in the interior, thus from a 
fraction to a large multiple of the possible life-time of the star. No complete 
mixing can be achieved by such a circulation. Not more than 0-2 of the inner 
mass can be assumed to get effectively mixed within a time interval of about 
3x 10° years, the rest being hopelessly separated; a stratification along the 
stream-lines themselves ensues, as the result of the changing composition and 
size of the core during such intervals of time (14). It may be inferred that in 
the polar direction heavier material containing less hydrogen will be sucked out 
from the core which gets more and more “exhausted’’; a systematic increase 
of the mean molecular weight from the equator toward the pole over the same 
level surface follows. This tends to diminish the oblateness in the central 
region and to increase it at the periphery*; with a sufficient contrast in molecular 

* This follows from 7 ~1/u over a level surface, whereas —d7/dn ~yg. 
20* 
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weight the variation of gravity over a level surface may become zero and, 
according to equation (13), the circulation may stop altogether at a certain distance 
outside the core; on the other hand, the rotational circulation will continue at 
an accelerated rate in the outer regions. ‘The conclusion is that the central 
region around the core can never be completely mixed with the outside by 
rotational circulation. 
"TABLE | 
Characteristics of Rotational Currents of Stellar Model ‘‘ Bj” 

Uniform composition with 0°03 metal, 0°57 hydrogen, o-40 helium content. 

M=2'117 x 10" gm.; R=6°30810!%cm. Convective core with o-1060 of total 

mass,7,;,=7°58xX10%cm. T,=2 107 deg.; log pp =2°1647; pe/pPym=72°6. Peripheral 

convective region from 7,=6:100 « 10!" cm. to surface, includes 1:13 x 10 -* of total 

mass. Lg=4°825 x 10° erg/sec. 

w=107* sec.~! 
-) wi los 
aa i eee r — log|te | log|tm| — (7) 
positive negative positive negative 

7°58 + 00 oo + 00 31-2 9°942 0°4780 

8'o 7'978 573 ° 32°8 8-438 04801 

8-8 7-508 587 0°2190 34°4 8-786 0°4844 

9°6 265 “128 0°3016 36°0 ‘O70 0°4913 
10°4 ‘100 073528 39°2 540 0°5146 
11°2 ‘979 *598 0'3886 42°4 956 0°5573 
12°0 ° 8-406 0°4138 45°6 0°6419 
12°8 8-847 8-275 0°4327 48°8 1'0448 
13°6 8-809 8-160 0°4476 49°32 - vee re) 
negative 
15°2 8-776 9°922 0°4676 52°0 6-725 
16°8 8-789 9°681 -0°4783 53°6 
18-4 8-829 9°424 0°4835 56°0 
20°'0 8-889 10°820 0°4852 58°4 
20°6 seo — 0°4853 59°2 

positive 

21°6 8-965 9°044 0°4850 59°6 
23°2 7048 9°468 0°4837 60-0 
24'8 7138 9°661 04817 60-4 
26°4 7°232 9°745 0°4796 60°8 
28-0 7*330 9°669 —0o'4780 61°00 
29°6 7426 T0°390 0°4774 
29°7 ove — 0 0°4774 


INN a] 
~ 
uv 
N 
nN 


x 

oo 
as 
wm 


x 


7 
i 
7: 
7° 
6: 
ra 


SIainints 


“il 


| ; vel wil wei 


+i | el nl 


Wi wi vol 





In the course of the circulation, especially at its initial stage (in which most 
of our chosen model remains even at the high assumed speed of rotation), the 
displacements are such that conservation of angular momentum leads to accelerated 
rotation in the equatorial plane; the amplitude in gravity is evidently increased by 
this effect, which thus tends to accelerate the circulation itself should the 
von Zeipel—Eddington theorem remain valid. However, as w is no longer 
constant, equation (1) cannot be strictly applied. In the time intervals considered 
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here, however, the contrast in molecular weight will be the more important 
factor and will stop the circulation in the interior in spite of the possible contrary 
effect of the equatorial acceleration of rotation. 

For a solar speed of rotation, w = 3 x 10-6, the time intervals given above are 
to be increased by a factor of about 1000; the circulation becomes so slow that 
its mixing effect during the life-time of a star is practically nil; the nuclear 
reactions in the core proceed much faster than the minute displacements produced 
by rotational circulation. The system of circulation differs in such a case from 
the preceding; no reversal of the direction of v, occurs because I > §a(p/p — 1) 
remains valid in (19) all through the radiative equilibrium zone. This zone 
constitutes one system of circulation, the central and the peripheral convective 
regions being linked together; however, the time of circulation from core to 
periphery, 3 x 10! years, is so large that the link has no practical significance. 

Taste II 
Stream-lines of Rotational Currents in Model “ Bj,”’ (14) 
w= 10~* sec. ~? 
r=radius in 10°cm.; ¢=latitude; t¢=time of complete circuit in years ; 
m, ~fraction of total mass of the model inside the stream-line surface of revolution 
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The internal law of rotation of the Sun is unknown; it is not impossible 
that w increases inwards but, as shown by the above computations, even a much 
larger value of w~10-* would not make much practical difference. If the 
suggestion, based on the evidence of lithium abundance, that the Sun to a 
considerable depth has gone through a phase of turbulent mixing at an early 
stage of its life corresponds to reality, it is very likely that uniform angular velocity 
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(as the result of eddy-viscosity) represents the initial state of solar rotation, except, 
perhaps, for a small metallic core of meteoric origin. 


The author acknowledges with gratitude the assistance given by Professor 'T’. G. 
Cowling in the form of critical comments and help in preparing the manuscript 


for press. 
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Summary 


A review of current theories of the scintillation of stars on optical wave- 
lengths shows that the refraction theory of scintillation requires excessive 
atmospheric density gradients, and fails to explain the observed scintillation 
in colour, whilst the physiological explanations are insufficient to explain all 
the observed phenomena. A new theory, in which the scintillation phenomena 
are attributed to Fresnel diffraction at a non-homogeneous atmospheric 
layer, is shown to explain the observed effects, and to require considerably 
smaller density gradients. A similar theory, taking into account the effect 
of a wave-length change of » 107, is shown to explain the observed fluctuations 
in intensity of radio waves from the localized extra-terrestrial sources of radio 


noise. 





1. Introduction.—Neither of the two main theories of the scintillation of stars 
on optical wave-lengths appears to be fully satisfactory. The physiological 
theory, as developed by Hartridge and Weale (1), fails to explain all the observed 
phenomena. These authors have shown that under certain circumstances the 


human eye registers scintillations even when illuminated by a constant-intensity 
source, and have suggested that the scintillation of stars may be due to this 
property of the eye. However, photographic observations by Gaviola (2), and 
photoelectric observations by Boutet (3), indicate that the intensity of the light 
from a star received by the eye does in fact vary, and hence that physiological 
factors cannot provide a complete explanation of the phenomena. 

The refraction theory of scintillation (2, 4) explains the production of light 
and dark striae some centimetres apart on the Earth’s surface by the focusing 
effect of irregularities of the same order of size and situated at a temperature 
inversion layer at a height of some 34km., but it will be shown that the theory 
requires excessively high density gradients. 

Fluctuations in the intensity of the radio emission from the localized sources 
of galactic noise were first observed by Hey, Parsons and Phillips (§) and recent 
work (6,7) has shown that these fluctuations are introduced by a local medium and 
may therefcre be analogous to scintillation on optical wave-lengths. 

In this paper the effect of Fresnel diffraction at a non-homogeneous trans- 
mitting layer is considered. It is shown that such a mechanism is capable of 
explaining the observed optical and radio phenomena by invoking known 
irregularities in the atmosphere and ionosphere respectively. 

2. The refraction theory of scintillation on optical wave-lengths.—T he refraction 
theory of the scintillation of stars on optical wave-lengths explains the production 
of striae on the Earth’s surface by means of the focusing effect of atmospheric 
irregularities at a temperature inversion layer (Fig. 1). 

* Originally submitted 1950 May 3. 
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An analysis of the theory, given below, shows that it may be criticized on three 
main accounts :— 

(a) It requires excessively high atmospheric density gradients. 

(6) The atmospheric ‘“‘lenses”’ postulated by the theory are so small that 
diffraction considerations show that they would have negligible focusing 
effect. 

(c) The theory fails to explain the observed scintillation in colour. 


The refraction theory as developed by various workers (2, 4, 10) postulates 
the presence of atmospheric “‘lenses’’, both convergent and divergent, spaced 
some 5 cm. apart at a range of about 3} km. from the observer. These will, of 
course, be irregular in shape, position and size, but as a simplification they may 
be idealized to a series of perfectly cylindrical lenses some 5 cm. apart, of radius 
of curvature ~ 2cm. and at a range of 3} km. 


| 
by 











Fic. 1.—The Refraction Theory of Scintillation (after Gaviola). 
XY represents the irregular boundary between two layers of air of different refractive indices, 
p, and ps. The boundary is taken to be at a height of approximately 34 km. and the irregularities 
in its surface, some 5 cm. apart, are thought to redistribute the light from a star by refraction into 
striae on the Earth’s surface. 


For such a lens we may write: 
I 


7a Ha) 5 


where f is the focal length of the lens, 7 is the radius of curvature (~2¢m.) and 
po and yu, are the refractive indices of the lens and the surrounding air respectively. 
In order that the focal length of the lens should be of the order of 3$km., it 


will be necessary that (4. — 44) = ort — p14) =3 x 10-8, since r =2cm. 


2 I 
r \3°5 x Io® 
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The refractive index of air at N.'T.P. is 1-0003 and changes by 1-0 x 10° 
for a 0-3 per cent change in density. Hence, if (4.—,,) is to be 3 x 10-8, a I per 
cent change in atmospheric density would be required in a distance of the order 
of 2cm. 

Such atmospheric density gradients are considerably in excess of those 
known to exist even under turbulent conditions at the Earth’s surface (8) and 
hence could not be expected to occur at heights of the order of 4km. It is of 
interest that this difficulty was first pointed out by Lord Rayleigh (9), whose 
comments, however, do not appear to have been considered in recent discussions 
on scintillation phenomena. 

A more fundamental difficulty is that diffraction considerations show that 
lenses of the size and range postulated could not focus the light effectively. ‘The 
resolving power of a cylindrical lens is given by A/a, where A is the wave-length 
and a is the lens aperture. Hence the image formed of a distant source will be 
of diameter D=2fA/a. Substituting the values of f, A and a postulated by the 
refraction theory, we have D = 10cm. for A =6000 A., and therefore the diffraction 
disk is larger than the lens itself, i.e. the lens is too small to focus the light 
appreciably. 

A third failure of the refraction theory is its inability to explain the observed 
scintillation in colour at high angles of elevation. ‘This has been attributed to 
the chromatic aberration of the atmospheric lenses (10). However, since 
(4 — 1) differs by only 1 per cent for red and green-blue light, the various colours 
will be concentrated to approximately the same degree by any given atmospheric 
lens, and scintillation in colour would therefore not be expected on a simple 
refraction theory. 

3. Some general considerations of diffraction at an irregular screen. (a) Intro- 
duction.—The discussion in the last section has shown that the atmospheric 
lenses postulated to explain optical scintillation phenomena are smaller than the 
first Fresnel zone in the disturbing region. ‘This fact suggests that the phenomena 
should best be considered as due to a Fresnel diffraction mechanism in which the 
idea of individual atmospheric lenses is replaced by the consideration of the 
distortion of the Fresnel Half Period Zone pattern of a plane wave viewed through 
a medium of slightly non-homogeneous refractive index. 

Such considerations indicate that the change in relative areas of odd and even 
numbered zones due to varying optical path lengths in the medium can markedly 
affect the observed intensity. Diffraction theory shows that the complete reversal 
of phase of one of the zones will increase the received intensity by 25 or 9 times, 
according to whether the zone is even or odd in number (Fig. 2() and (c)). 
Such a phase change could hardly occur from a random distribution of refractive 
index irregularities, but a change in the relative area of the odd and even numbered 
zones, at any particular instant, of about 0-2 of a zone would not be unlikely. 
In such a case the received intensity would be increased or decreased by about 
three times, depending upon the area affected and the change in phase (Fig. 2 (d) 
and (e)). 

It is not practicable to pursue these simple Fresnel zone considerations to 
explain the scintillation phenomena due to an irregular diffracting screen, since 
it would be necessary to sum the effect of each of the resultant distortions of the 
Fresnel zone pattern over the whole wave-front. In this case statistical methods 
must be used, such as those developed by Booker, Ratcliffe and Shinn (11). 
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These authors have shown that the effect of an irregular diffracting screen may be 
analysed in terms of the generalized auto-correlation function of the electro- 
magnetic field strength across a plane immediately below the disturbing layer. 
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(a) Phase change=o. (6) First zone changed by 7. (c) Second zone changed by 7. 
No resultant change in amplitude. Amplitude increased three Amplitude increased to 
times, intensity nine times. five times, intensity 25 

times. 





A 


(d) Central portion of first zone advanced by  (e) Central portion of first zone retarded by 
¢<n/2. ¢ < 2/2. 
Amplitude increased to 1-7 times, intensity to Amplitude decreased to 0-6 times, intensity to 
approximately three times. 0°36 times. 
Fic. 2.—Effect of Distorting the Fresnel Zone Pattern. Amplitude-Phase Diagrams. 

AB= resultant of first zone. 

BD=resultant of second zone. 

DE =resultant of all subsequent zones (assumed undistorted). 

AE=resultent of undistorted pattern. 
AE'=resultant of distorted pattern. 


It has been proved (11) that the average angular power spectrum observed 
by a recording device receiving radiation from an irregular diffracting screen is 
the Fourier transform of the average auto-correlation function p,(&)* of the 
wave-field across a plane immediately below the disturbing region. This means 
that an observer receiving radiation through an optical screen possessing irregular 
transmission properties receives radiation not only along the line of sight but 
from a cone whose angle is determined by the auto-correlation function p,(&), 
i.e. by the ‘‘ average size’’ of the irregularities in the wave-front. 


* See appendix for definition of the auto-correlation function. 
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For the optical and radio wave-lengths considered in the present paper, 
absorption or reflection in the atmosphere or ionosphere can be neglected, and 
the amplitude of the wave-field immediately below the disturbing layer may be 
taken as constant. However, due to the differing path lengths of the radiation 
through the disturbing regions, the emergent wave-front will be distorted and we 
now proceed to discuss the effects of this distortion at different distances below 
the disturbing layer. In this connection it will be found that the ratio of the radius 
of the first Fresnel zone, 1/(RA), to the distance £ at which p,(€) falls to zero is 
an important criterion determining the type of scintillation phenomena observed. 

(6) ./(RA) large compared with £.—It is proved in the appendix that when 
(RA) = &, the angle subtended by one irregularity in the emergent wave-front 
is approximately equal in size to the received cone of radiation. By the property 
of Fourier transforms, if £ is now taken to be small compared with +/(RA), the 
received cone of radiation will be large compared with the angle subtended by 
the first Fresnel zone. But since the irregularities in the emergent wave-front 
are small compared with the first Fresnel zone, this means that a large number* 
of irregularities are contributing to the observed intensity at any given point. 
These irregularities have random phase and hence the resultant intensity will have 
a Rayleigh distribution. 

A further point of interest in this case is the effect of varying the observing 
wave-length slightly, on the assumption that the phase change introduced by an 
irregularity in the disturbing region is independent of wave-length. The phase 
of a point on the wave-front at the extreme edge of the received cone relative to 
the centre point is critically dependent upon the wave-length (owing to the large 
width of the cone relative to the first Fresnel zone) and it can be shown, on the 
assumption that the irregularities are large in size compared with A, that if N is 
the ratio of the radius of the first Fresnel zone to the size of the irregularities, 
a change in the wave-length of A/N® is sufficient to produce an uncorrelated 
scintillation record. 

In the optical case, where the human observer uses a band-width of approxi- 
mately one octave, this difference between the scintillations on different wave- 
lengths could be an important factor limiting the degree of scintillation observed. 
For instance, if the scintillations cease to correlate for differences in frequency of 
only 0-OI octaves, approximately 100 independent scintillation records would 
be observed simultaneously and the total resultant degree of scintillation would 
therefore be reduced to approximately one-tenth of that on a single wave-length. 

(c) (RA) small compared with &.—In this case the Fourier transform of 
pz(€) is small compared with that of the first Fresnel zone, and hence the received 
cone incident at a point on the Earth’s surface is smaller than that subtended by 
the first Fresnel zone and therefore is considerably smaller than that subtended 
by the irregularity itself. ‘This means that the wave-front is almost plane over 
the acceptance cone and hence the degree of scintillation will be very small. 
Also, owing to the very small variation in phase across the cone there will be a 
considerable degree of correlation between the scintillations on different wave- 
lengths, and therefore scintillation in colour would not be observed. 


* Approximately N*, where N= +/(RA)/éo.. The angular width of the received cone will be N 
times larger than that of the first zone and therefore covers N* x the area. But there are approxi- 
mately N® irregularities per zone and hence the total number of effective irregularities is 
approximately N‘. 
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(a) Diffraction from an irregular screen. 


The width of the scattered cone of radiation is 
determined by the auto-correlation function 
p,A€) across the emergent wave-front. 
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(c) V/ (RA) < £o. 


The scattered cone of radiation is small compared 
with the angle subtended by one irregularity 


in the emergent wave-front. 


(b) +/(RA) > £0. 
The scattered cone of radiation is large compared 
with the angle subtended by one irregularity 

in the emergent wave-front, 
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(d) 1/(RA)=£o- 
The scattered cone of radiation is approximately 
the same size as the angle subtended by one 
irregularity. 


Fic. 3.—Diffraction from an Irregular Screen. 


N.B.—The distortion of the emergent wave-front has been considerably exaggerated. 
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(d) /(RA) approximately equal to €,.—In this case, as shown in the appendix, 
the received cone is approximately that subtended by an individual irregularity. 
The intensity observed will therefore depend critically upon the nature of the 
curvature of the wave-front inside the acceptance cone. ‘The degree of scintilla- 
tion will be considerable, but owing to the narrowness of the scattering cone, 
there will be a higher degree of correlation between the scintillation phenomena 
on different wave-lengths than in the case for ,/(RA) large compared with &,. 
The limit in the optical case will be imposed by the fact that the phase change 
introduced by the irregularities is not constant with wave-length but is propor- 
tional to 1/A, and also that the width of the acceptance cone is proportional to A. 

(e) The dimensions of the Fresnel diffraction pattern.—It is also proved by 
Booker, Ratcliffe and Shinn (11) that the generalized auto-correlation function 
across the Fresnel diffraction pattern (i. e. across the scintillation pattern observed 
at ground level) is the same as that across the emergent wave-front itself. ‘This 
result means that by observing the dimension of the diffraction pattern at ground 
level we are able to deduce the size of the distortions in the wave-front as it emerges 
from the disturbing region. 

The observation that the “‘striations’’ in the scintillation patterns are some 
5 cm. apart means therefore that the phase of two points more than 5 cm. apart on 
a plane immediately below the disturbing region is not correlated, i.e. that the 
distance over which the emergent wave-front is plane (to better than A/4) is of 
the order of 5 cm. 

4. The application of the diffraction considerations to optical scintillation. 
(a) The range of the disturbing region.—The observations of Respighi (12) and 
Ellison (10) on colour scintillation indicate that there is in general no close corre- 
lation between the scintillation patterns of light of wave-lengths differing by more 
than about 1000Aa. ‘This fact therefore suggests (by the reasoning of Section 3 
(a), (6), (c), (d)) that the irregularities in the emergent wave-front are of the order 
of size of the first Fresnel zone \/(RA). But by the theorem of Booker, Ratcliffe 
and Shinn quoted in Section 3(e), the observation that the “striations’’ in the 
scintillation pattern observed at ground level are some 5cm. apart must mean 
that the irregularities in the emergent wave-front are some 5 cm. in size. Hence 
by writing \/(RA)>5cm. we find that the disturbing region is at a distance of 
approximately 4km., and that the angle subtended by these irregularities in the 
emergent wave-front is approximately 2} seconds. 

(6) Physical requirements of the diffraction theory.—T he total increase in optical 
path length introduced by the presence of the Earth’s atmosphere is about 
4 x 108 A at normal incidence for A=6000A. The above result therefore indicates 
that the path change introduced by the atmosphere varies irregularly by about 
one part in 107 for two points 5 cm. apart at the Earth’s surface. ‘These variations 
may have been introduced by the passage of the light through an individual, 
small, relatively strong irregularity, or through a number of weaker irregularities. 
If the irregularities are assumed to be of the order of 10 cm. in size (8) the density 
gradient required is of the order of 0-003 per cent in 5 cm. for a disturbing region 
o-2km. thick. Such gradients (approximately one-thousandth of those required 
for the refraction theory) are of the same order of magnitude as the mean vertical 
density gradient in the disturbing region, and are about one-fiftieth of those 
observed under turbulent conditions at ground level. 

(c) The duration of individual scintillations.—As the star moves across the sky, 
the cone of light received at a single point on the ground moves relative to the 
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optical depth irregularities. The ratio of the angular subtension of the 
irregularities in the wave-front to the angular velocity of the star gives an estimate 
of the duration of the scintillations. For stars of declination 0° the rate of motion 
is a maximum and is about 15 seconds per second. Since the angular subtension 
of the irregularities is about 3 seconds (Section 4(a)) this suggests that the 
scintillation period should be of the order of 0-2 second. However, winds 
having a component of velocity 30 km./hr. in the plane of the emergent wave-front 
will reduce the scintillation period to about 0-01 second. Owing to the persis- 
tence-of-vision effect in the retina of the eye only the frequency components of 
the scintillations of less than approximately 30 cycles per second will be observed 
by the eye. This may well be the explanation of the comparative absence of 
scintillation on a windy night, although it should be borne in mind that a strong 
wind might tend to destroy the irregularities responsible for scintillations. 

(d) Scintillation in colour.—'The diffraction considerations described in 
Section 3 introduce colour effects in two ways :— 


(i) The angular width of the cone of light incident at a given point on the 
Earth’s surface is proportional to the wave-length. 


(ii) The phase change introduced by a given irregularity is proportional to 
1/A. 


These two effects combine to cause light of different wave-lengths to be 
distorted in different ways, and hence to make the star scintillate in colour. 
This is true of all angles of elevation. In addition, however, there is a third 
effect, important only at low angles of elevation. For stars at elevation lower 
than about 15° the incident cones of radiation for red and blue light are separated 
by more than 2 seconds (owing to dispersion effects in the atmosphere) and hence 
the two cones are affected by different atmospheric irregularities. For this 
reason we would expect scintillation in colour to be more marked when the 
sources are at low angles of elevation. 

(e) Variation of scintillation activity with angular subtension of source.—If the 
source subtends an angle very much less than that subtended by one irregularity, 
each point of the source will be affected similarly by the irregularity and the source 
will therefore be observed to scintillate as a whole. If, however, the source 
subtends an angle large compared with that subtended by an irregularity, then 
the different parts of the sources no longer scintillate in the same manner, an 
increase in the intensity from one point tends to be compensated by a decrease 
from some other point, and the source will therefore not be observed to scintillate. 

This result is well borne out by the behaviour of the planets whose semi- 
diameters range from about 1 second to about 30 seconds. The smallest ones 
are not visible to the human eye and scintillation information is scanty, but those 
of middle angular subtensions are rarely found to scintillate and the larger ones 
very infrequently or never. 

(f) Changes in the apparent position of the source.—Due to the phase- 
distorting irregularities in the disturbing region, the phase-gradient across an 
observing plane will not be uniform and the apparent position of the star will 
vary irregularly. ‘The amplitude of these angular deviations will be given approxi- 
mately by (A/2)/f, i.e. by the ratio of the deviation of the wave-front from a 
plane to the linear dimensions of the irregularities. On optical wave-lengths 
the apparent position of the star would therefore be expected to change 
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irregularly, at the scintillation frequency, by ~1 second of arc. This is in 
good agreement with experimental observations (4). 

(g) The production of striations.—As described in Section 3 (e) the theoretical 
work of Booker, Ratcliffe and Shinn shows that the generalized auto-correlation 
function of the Fresnel diffraction pattern is the same as that of the emergent 
wave-front. If now for some reason (e.g. an atmospheric wind) the atmospheric 
irregularities are not isotropic, the generalized auto-correlation function of the 
diffraction pattern at ground level will be a function of direction across the ground. 
In cases where the auto-correlation function varies considerably with direction 
the diffraction pattern will take the form of an approximately parallel series of 
bands or “‘striae’’ of light and shade, instead of an irregular series of randomly 
shaped and orientated patches. 

5. Comparison of the refraction and diffraction theories of the scintillation of 
stars on optical wave-lengths.—T he two theories may be summarized and compared 
by means of the two tables given below. ‘Table I summarizes the requirements 
of the two theories, whilst Table I] summarizes the manner in which they explain 
the observed phenomena. 


TABLE | 
Physical Requirements of the Refraction and the Diffraction Theories 

Assumptions of Theory Refraction Theory Diffraction Theory 
Range of disturbing region Approx. 3°5 km. Approx. 3°5 km. 
Atmospheric density gradient 0°5 per cent per cm. 0'0006 per cent per cm. 
Equivalent atmospheric temperature 1°5 deg. C. per cm. 0002 deg. C. per cm. 

gradient 
Size of irregularities Approx. 5 cm. Approx. 5 cm. 


TABLE [1 
Explanation of Ob »ved Phenomena on the Refraction and the Diffraction Theories 


Explanation on the Explanation on the 


Observed Phenomena ‘ Saige > ; 
Refraction Theory Diffraction Theory 


Scintillation in colour Not explicable, except for Explains _ scintillation 
stars very low on the colour at all angles 
horizon. elevation. 


Variation of degree of Explains these pheno- Explains scintillation pheno- 
scintillation with size of mena as due to re- mena as due to diffraction 
source fraction by atmospheric by atmospheric irregular- 
waves at the narrow ities in a_ relatively 
Duration of individual boundary of two atmo- thick disturbing region. 
scintillation spheric layers of different Requires normal atmo- 
density. Requires ab- spheric density gradients. 
Changes in apparent posi- normal atmospheric 
tion of the stars density gradients. 


6. The application of the diffraction considerations to radio wave-lengths. 
(a) Introduction: the experimental observations of scintillation effects on radio 
wave-lengths.—F luctuations in the received intensity of the radio emissions of 
extra-terrestrial origin were first observed by Hey, Parsons and Phillips (5). 
The subsequent discovery of localized sources of radio emission (or radio stars) 
(16, 17) led to the general belief that such fluctuations were due to the variations in 
the emission from the sources. Recent work by Smith (6) and Little and Lovell (7) 
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has shown that receivers more than 20 km. apart on the Earth’s surface register 
different scintillations and hence the fluctuations must be locally impressed. 
The suggestion that these fluctuations could be explained in terms of the dis- 
tortion of the Fresnel zone patterns by irregularities in the ionosphere was made 
by Little and Lovell (7) and the application of the diffraction theory is here given 
in more detail. 

(b) The range of the disturbing layer.—Experiments, to determine the correla- 
tion between radio scintillations observed at one site on different frequencies, have 
shown that in general there is no correlation between records taken at frequencies 
differing by 30 per cent (13). More recent work (20) has indicated that good 
correlation exists up to frequency differences of the order of 5 per cent to 10 per 
cent. In view of the diffraction considerations of Section 3, these results suggest 
that the irregularities in the emergent wave-front have dimensions of the order 
of the first Fresnel zone at the range of the emergent wave-front. Spaced 
receiver observations (6, 7) have shown that the correlation between fluctuation 
records at the same frequency remains good up to spacings of about 4km. and, 
by the second theorem of Booker, Ratcliffe and Shinn (Section 3 (e)), this means 
that the irregularities in the emergent wave-front are approximately 5 km. in size. 
Hence, equating 1/(RA) to 5km., we obtain a figure of R= 4o00km. for A=6 
metres. However, this relationship between \ (RA) and &, is only an approximate 
one and, owing to the square root term, the figure for the range of the disturbing 
region is liable to considerable inaccuracy. In order to determine the probable 
location of the disturbing region more exactly, it is necessary to make use of other 
information, such as that obtained by correlating the occurrence of galactic noise 
fluctuations with anomalies in the various ionospheric layers. 

(c) Correlation of scintillation activity with the presence of F. ionospheric 
srregularities.—An analysis of the occurrence of galactic noise fluctuation at the 
transit of the Cygnus and Cassiopeia point sources, during the period 1949 May 1 
to 1950 November 30, and of the occurrence of ‘‘spread F”’ conditions in the 
F, layer is given elsewhere (14). It is sufficient to state here that a mutual coin- 
cidence of at least 75 per cent exists between the occurrence of radio scintillations 
(at transit) and the observation (by ionospheric workers) of ‘‘ spread” or “‘ diffuse” 
echoes from the F, layer. ‘These echoes are thought to be due to the presence 
of electron clouds of varying density and height in the layer. (No similar 
significant correlation has been found with sporadic E phenomena.) 

This high degree of correlation suggests that the galactic noise fluctuations du 
indeed originate in the F, layer. If we substitute the height of the maximum 
electron density in the F, layer (approx. 400 km.) we find, using a figure of 5 km. 
for €, that 1/(RA)/& varies between 0-3 for a source in the zenith to 0-7 for a 
source at low angles of elevation. ‘The diffraction considerations of Section 3 
indicate that these values for the ratio \/(RA)/£ are sufficiently near to unity 
to explain scintillation phenomena of the type observed. 

(d) Physical requirements of the diffraction theory to explain the observed radio 
scintillation phenomena.—The spaced receiver observations already referred to 
suggest that the ionospheric irregularities are of the order of 5km. in size. Since 
scintillation phenomena will not be observed unless the irregularities possess 
sufficient refractive index changes to increase or decrease the radio path length 
by about half a wave-length, it is possible to estimate the minimum change in 
electron density required to produce scintillations. 
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The refractive index of an ionospheric layer to a radio wave of frequency 
f is given by 
. Ne* 
lifes _— amf? ’ 
where N=number of electrons/c.c., m=mass of electron, and e=charge of 
electron. 


If f, is the critical frequency corresponding to an electron density N, we have 


af Gy 


In the case under discussion f is large compared with f,, and hence we may 


write 
pat~ i(f./f)?- 
If Au is the difference in refractive index between an ionospheric irregularity 
and the surrounding ionosphere, 
wee Rm 2_ f2 
Apu pi aplli fo); 
where f, is the critical frequency of the irregularity itself. ‘The change in radio 
path length, introduced by an irregularity of linear dimension L and refractive 
index differing by Ay from the surrounding medium, will be LAyw/A wave-lengths. 


If this is to equal }A (the requirement for the production of scintillations) we 
have 


L\ 1 p 
eh Woe 2_ #2) 1 
(x) 2f* (fi fo) - 
‘Taking A=6 metres, L =5000 metres and f, = 3°0 Mc./s., we have 


par 


fi= Fz t+fe=12 x10" or f,=3'5 Me./s. 


If Ap is taken as negative, we obtain f, 2-5 Mc./s. ‘The corresponding electron 
densities are 7:5 x 10‘ electrons/c.c. and 1-5 x 10° electrons/c.c. ‘These figures 
are based on the assumption that the disturbing region is only one irregularity 
thick. If there is more than one irregularity in the line of sight, smaller changes 
in electron density would be sufficient to produce the required distortion of the 
emergent wave-front. Investigations of the F region of the ionosphere indicate 
that under strong spread F conditions the critical frequency can extend over a 
range of two or more Mc./s. (1g). Such a range would appear to be more than 
sufficient to produce radio scintillations. 

(e) The durations of individual fluctuations.—As in the optical case, it is possible 
to explain the observed durations in terms of the ratio of the angular subtension 
of the irregularities to the angular velocity of the source; and of the ratio of the 
size of the irregularities to the velocity of the wind in the disturbing region. 
The angular subtension of one irregularity in the emergent wave-front varies 
between some 40 minutes (in the zenith) to about 10 minutes (at low angles of 
elevation). Since the angular velocity of the source is of the order of 10 minutes 
of arc/minute, fluctuations of period about 4 minutes would be expected at 
transit, in the absence of any appreciable ionospheric wind effects. However, 
winds of the order of 1000 km./hr. in the F, region of the ionosphere have been 


21 
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observed (15) and would therefore produce fluctuations of period approximately 
20 seconds. These figures agree well with the observed range in durations of 
from approximately 10 seconds to about 5 minutes for fluctuations at transit (14). 

(f) Variation of scintillation activity with angular subtension of source.— 
Following a similar argument to that used for the optical case, sources whose 
angular subtensions are large compared with 40 minutes of arc will not be observed 
to scintillate. The radio stars have been shown by Bolton (16) and by Ryle 
and Smith (17) to be smaller than about 6 minutes of arc and hence all these sources 
would be expected to scintillate. Owing to the great difficulty of observing all 
but the four most intense radio stars, scintillation information on other radio stars 
is scanty, but there is as yet no well authenticated case of a source which does 
not scintillate. 

(g) Variation in the apparent position of the source.—As in the optical case, 
the apparent position of the star will tend to fluctuate at the scintillation frequency 
through an angle of about one-third the width of the scattered cone. For 5km. 
irregularities and a wave-length of 4 metres the scattered polar diagram will have 
a first minimum at an angle of approximately 3 minutes of arc from the direction 
of the incident radiation. In this case the random fluctuations in the apparent 
position of the source would be expected to have an amplitude of the order of 
I minute of arc. ‘This figure is in good agreement with the observations of 
Ryle and Hewish (18) who have shown by using a high-resolution radio-inter- 
ferometer technique that the random changes in apparent position are of the order 
of 1-2 minutes of arc. 

(h) Variation of degree of scintillation with wave-length.—The phase change 
produced by a given irregularity is proportional to the wave-length and hence 
the degree of scintillation would be expected to decrease with decreasing wave- 
length. This effect has been observed (13), scintillations being more intense 
on the longer wave-lengths. It should however be pointed out that any theory 
seeking to explain the scintillations in terms of refractive index irregularities in 
the ionosphere would predict the decrease of scintillation with decreasing wave- 
length, since for the ionosphere (4 —1) is proportional to A*. 

7. Conclusion.—A new theory of scintillation, based on Fresnel diffraction 
at a non-homogeneous layer, has been shown to be capable of explaining the 
observed phenomena on optical and radio wave-lengths. In the optical case 
the theory requires considerably smaller atmospheric density gradients than are 
required by the simple refraction theory. ‘The use of photoelectric cells and a 
photographic recording system should enable the application of the diffraction 
theory to optical wave-lengths to be more rigorously checked. On radio wave- 
lengths the study of the scintillations of the radio stars should provide a new method 
for studying the complexities of the F region. 

8. Acknowledgments.—The author wishes to express his sincere thanks to 
Dr A. C. B. Lovell, Director of the Jodrell Bank Experimental Station, for his 
constant help and encouragement during the preparation of this paper; to 
Mr J. A. Ratcliffe of the Cavendish Laboratories, for invaluable criticisms and a 
most helpful discussion of the application of statistical methods to diffraction 
at an irregular screen; to Dr M. A. Ellison, for his helpful criticism and interest 
in this work; and to Mr A. Maxwell, for assistance in taking the radio obser- 
vations. He is also indebted to the Department of Scientific and Industrial 
Research for the award of a maintenance grant. 








No. 3, 1951 _ scintillation of stars on optical and radio wave-lengths 


APPENDIX 
The auto-correlation function and its Fourier transform 
‘The purpose of this appendix is: 
(a) To define the auto-correlation function p,(é), and 
(6) To show that the angle subtended by one irregularity in the emergent 
wave-front is approximately the same size as the scattered cone of 


radiation when &, = \/(RA) (i.e. when the size of the distortions in the 
emergent wave-front is about that of the radius of the first Fresnel zone). 


The terminology throughout is that used by Booker, Ratcliffe and Shinn (11). 


t. The auto-correlation function.—Consider a plane wave incident upon an 
irregular diffraction screen which varies in one direction only, the X direction. 

Let E,(x) be the function describing the amplitude distribution, in the 
direction X, of the X component of the electric field across a plane immediately 
below the disturbing region. ‘Then the auto-correlation function, in the direction 
X, of this component of the electric field is given by 


[°° E&(x). B(x +8).dx 
| j : E* (x). E,(x). dx : 


. 





pil€) = 


Here E*(x) is the complex conjugate of the amplitude of the Y component of the 
electric field at the point x and E,(x+£) is the (complex) amplitude of the X 
component of the electric field at a point (x +), all distances being measured in 
wave-lengths. ‘This is a symmetrical function, having a maximum value of 
unity for =o, and falling off towards zero as € increases. &,, the distance at 
which the auto-correlation function has fallen to zero, is a measure of the average 
size of the irregularities in the emergent wave-front. 

2. The ratio of the angle subtended by one irregularity in the emergent wave -front 
to the semi-angle of the scattered cone of radiation.—\t is shown by Booker, Ratcliffe 
and Shinn (11) that the average angular power spectrum in the scattered cone of 
radiation is proportional to the Fourier transform of the auto-correlation function 
defined as above. If now we assume that this auto-correlation function is of 
the Gaussian error function type, the Fourier transform of the auto-correlation 
function will also be Gaussian, and we have 


prl(€) =e sd aallt | P(s) Pae27*a"s" 


where a is the average number of wave-lengths separating two points in a plane 
below the disturbing region between which the correlation is 0-61, and 1/(27a) 
is the sine of the angle at which the scattered radiation has fallen to 0-61 in power. 

Multiplying 1/(27a) by the range of the disturbing region R and equating 
to ad (the condition for the average irregularity to subtend an angle equal to the 
scattered cone of radiation) we have 


ramon a-(2)V(8) 


If we take &, the distance at which the auto-correlation function has fallen to 
zero, as 2aA, then 
9 =0°84/(RA). 
21* 
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Hence we see that the scattered cone of radiation is approximately equal to the 
angle subtended by one irregularity in the wave-front when = 4/(RA). When 
£ is small compared with 1/(RA), the scattered cone of radiation will be large 
compared with the angle subtended by 1/(RA) and will therefore be very much 
larger than the angle subtended by one irregularity; on the other hand, if &, is 
large compared with 1/(RA), the scattered cone of radiation will be very small 
compared with the angle subtended by one irregularity. 


Jodrell Bank Experimental Station, 
University of Manchester : 
1950 December 21. 
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Summary 


A theoretical study is made of the rate at which free atoms form molecules 
by radiative association. Previous treatments of this problem are shown to 
be in error. A revised formula for the coefficient expressing the rate is 
derived. Approximate methods of evaluating this formula are discussed; 
and tables, designed to ease the computational labour, are presented. 

The methods developed are used to obtain the rate coefficients for CH, 
CH~, Nf and Hf formation as these are of astrophysical or geophysical 
importance. In the case of Hj the calculations can be performed with some 
precision. Free-free radiative collisions between H atoms and H* ions are 
also investigated. 





1. Introduction.—At low gas densities molecules are formed mainly by 
radiative association 


X+Y—>XY+hy, (1) 
since there are too few third bodies for the reaction 
X+V+Z—-XV+Z (2) 


to compete. ‘lhe density below which (1) is more rapid than (2) is dependent 
on the atoms involved, and no precise general statement regarding it can be 
made. But unless the transition which stabilizes the system is optically forbidden, 
the critical concentration of Z is unlikely to be less than 10"/cm.®, and is frequently 
several orders of magnitude greater. Consequently radiative association must be 
the dominant process under many naturally occurring conditions. The rate at 
which it proceeds can be expressed in terms of a coefficient y, defined by the 
equation 


dn(X Y )/dt = yn(X)n(Y), (3) 


where the n’s denote the concentrations of the species indicated. In an important 
recent paper Kramers and ter Haar (1) studied the problem of determining the 
coefficient. Unfortunately an error occurs in the classical part of their treatment, 
as a result of which the influence that a short-range attractive interatomic field 
has in increasing the number of collisions was seriously overestimated; and in 
the quantal part the possible implications of certain correlation rules were not 
fully realized. Radiative association is re-investigated in the present note. 
A general formula for its coefficient is derived. This is then applied to compute 
the rates of formation of CH,CH*+, N} and Hj} (which are of interest to astro- 
physicists and geophysicists). 
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2. General formula 

2.1. Suppose that two normal atoms, whose internuclear distance is r, 
approach each other. ‘The fraction of occasions on which they move along a 
particular energy curve, U,(r), may generally * be taken to be g, the ratio of 
the statistical weight of this curve to the sum of the statistical weights of all 
possible curves. If the Einstein transition probability to a lower potential 
energy curve U,(r) is A(r), then the chance of a molecule being formed is 


g | A(r) dt, (4)t 
the integration being carried out for all time ¢ for which r is such that the emission 
of a photon deprives the system of so much energy that the atoms cannot separate 
indefinitely. ‘To obtain the rate coefficient associated with the process it is 
necessary to average over all values of the relative velocity v and of the impact 
parameter p. By the laws of motion dt and dr are connected by the relation 


pv? = 2U;,(r))\-1? 
sti Ee a ee 
dt =dr {° x - , (5) 


where m is the reduced mass. Using this and assuming a Maxwellian distribution 
at temperature 7’, it may be shown that 


32 r¢ 2 2ey2 y —1/2 
y =8n'*g (sir) | | | [ exp (- =) vA(r)p {e - a - aad ] dv dp dr. 
(6) 
k being Boltzmann’s constant (1). Physically the most natural procedure is to 
calculate the contribution from each orbit, that is to integrate over the range 


of r appropriate to each p and wv value. But as Kramers and ter Haar (1) have 
pointed out, it is mathematically easier to carry out the p integration first. 
Obviously the lower limit is simply zero. ‘The upper one is py, the greatest 
value of p which is consistent with the chosen v and r._ Integration yields 


ft ae mv" 
y= Sng (—F 7) | | Aertvexp (- =r) 


x {(v? —2U(r)/m)!? — (v* —2U,(r)/m — piye*7*)} dr de. (7) 


At a perihelion point dr/dt vanishes. From (5), therefore, earlier workers have 
supposed that p,, is given by the equation 

Pin = {7? -2U,(r)r?/mo*}, (8) 
and hence that the second term in (7) is zero. This supposition is in general 
incorrect. Consider the graph of {r?—2U,(r)r?/mv*} versus r*. If U,(r) falls 
off more slowly than as r* the resulting curve increases monotonically; no 
complexity arises and (8) is in fact valid. But if (as is commonly the case) U,(r) 
falls off more quickly than as r~? the situation is entirely different. ‘The curve 
passes through an initial maximum, then through a minimum, and finally it too 
increases monotonically (cf. Fig. 1). In the regions r<(OA)? and r>(OC)!” 
(8) can again be applied. However, in the important intermediate region, 
(OA)! <r<(OC)#", naive use of the equation is erroneous. Thus when r is, 
for example, (OB)1?, it gives p,, to be (BE)"*; and if p had this value the atoms 
could not reach a distance (OB)? from one another, since dr/dt would vanish at 


* See however Section 3.1. 


+ It is assumed here, and throughout the paper, that the radiation field is so dilute that 
stimulated transitions are unimportant. 
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a separation of (OD)#*. Clearly py, is actually only (CF), as this is the greatest 
value of p for which the approach is possible. 

The variation of U,(r) is usually extremely rapid. Consequently, unless T 
is very high, it is permissible as a crude first approximation to take (CF)"? to be 
independent of v and equal to an effective collision radius p, such that 


Si U,(p) rig tkT, (9) 


that is such that the potential energy at this internuclear distance is half the mean 
kinetic energy of approach. 








0 iA B ic D 











Fic. 1.—Schematic graph of {r?—2U,(r)r*/mv*} (vertical scale) versus r* (horizontal scale). 


For the integration over v the lower limit v, is (2U,(r)/m)"*, where 
U,(r) is positive, and zero where U,(r) is negative*; the upper limit v, is 
1 P v7 es Pp 2 

{2(U,(r) — U,(r))/m}?. Elementary analysis yields 


y=4ng(P— 7), (10) 
where P= ¥ A(r)r?{ — F (x1) + F(x»)} exp (x7) dr, (11) 


“¢ 


| AG yP(Hxs)— Ha} (6 
] 
| 
J 


2 
> 


:. exp (— x3) dr, (12) 


: | xtexp( — x") dx, (13) 
vB (x real), 
(14) 
(x imaginary), (15) 
x? = — U,(r)/kT, x3= — U,(r)/kT, 


ti=at/(G-1), st-s}—sde at (16) 


[= exp (x) dx, 





* It is assumed here that at small r, U,(r) is positive, and that at moderate and large r, it is 
negative (tending towards zero at infinity). Other cases are discussed in Section 2.4. 
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8’ being the value of r at which U,(r) changes sign, 5” being that at the inner 
boundary of the intermediate region where (8) is inapplicable ((OA)** of Fig. 1) 
and p being the effective collision radius already defined. The r integration 
(with the limits as indicated) covers only half the path; to compensate for this 
an external multiplying factor of 2 has been introduced. Kramers and ter Haar 
derived the same expression for the Y term but did not take into account the 
@ term (which may be of almost equal magnitude). 

The difference between 5’ and 8” is negligible in most instances so that the 
distinguishing suffixes may be omitted; except at very high temperatures it is 
usual that — U,(r), — U,(r) and U,(r)— U,(r), >kT for the more important part 
of the region 5<r<p; further, the condition — U,(r)>kT generally remains 
valid for that part of the region beyond in which A(r) is appreciable. In these 
circumstances considerable simplification can be effected. ‘Thus on making use 
of the asymptotic expansions for the F and Y functions, (10) reduces to 


(kT \'? /° ‘ ere 
y=ane| p*(—) | Ary —Uyy drs |” A(rye— F(x) exp (a4 dr |. (27) 
Standard tables can be employed in the evaluation of A(x,), since of course 


F (x) =9(x) - BA xexp (— x), (18) 


il 


where (x) is the error integral. In fact, however, serious loss of accuracy is not 
caused by replacing the second term in the brackets by 


i: A(r)r* dr. (19) 


2.2. ‘The development just given is useful in that it shows in a simple fashion 
the type of correction that must be applied to the original formula of Kramers 
and ter Haar. But the approximation that p,, is a constant is unsatisfactory 
and unnecessary. At large internuclear separation it is usually possible to 
express the potential by an inverse power law 


U,(r) =ar-’. (20) 


If this is permissible for part of the region where — U,(r)>T, then for all values 
of v of significance the minimum value of the function r°(1—2U,(r)/mv*) is of 
magnitude 


{A/(A—2)}{ — a(A—2)/mv?}*, (21) 
and occurs at an internuclear distance of 
{ —a(A—2)/mo*}"". (22) 
Thus a given r lies without the position of the minimum if 
v? > —ar-(A—2)/m, (23) 
within if 
uv? << —ar~“(A—2)/m, (24) 
and is less than 5” (defined as above) if 
r(x —2U,(r)/mo*) <{A(A—2)}{ — a(A— 2) /mo®}*”. (25) 
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Consequently pj, is given by (8) if (23) or (25) is satisfied, and by (21) if (24) is 
satisfied. Hence (7) becomes 


y =8rltg | i JI (r)A(r)r? dr, (26) 


where 


I(r)= | {(e+H)®—U}exp (0) de, 


U =o (for the ranges of v given by conditions (23) or (25)) 


=e+H- = {(A —2)K/2e}*” (otherwise), 
€=mv*/2kT, 
H = —U(r)kT, 
K=—ar ‘kT, 


and where the limits «, and «, correspond to the v, and v, of the other approximation. 

As mentioned before, 5’ and 5” are usually almost equal, so that the region 
betweenthem is unimportant. There would of course be no difficulty in evaluating 
the contribution from it should the requirement arise. The integration over the 
remaining range is simplified by considering two extreme cases. 

Case 1.—lIf H>1, U is given by (29) except for velocities well beyond the 
maximum of the Maxwellian distribution. ‘The difference between the two 
square roots in (27) can be obtained by the binomial expansion. The first 
non-vanishing term yields .4(r)=.4(r)*, where 


I id 2)) ‘te jaya 

JI(r)* = se {A/(A—2)}{ K(A — 2)/2}" 3 «exp (—e) de. (33) 
Since U,(r) is negative, «, is zero. If in addition, U(r) — U,(r) >RT, €, may be 
taken as infinite, so that 

I 


I(r)4 = a5 {K(A—2)/2}"P((A—2)/). (34) 


2H 
The retention of the next term of the expansion gives, as an alternative approxi- 
mation, 


I(r)! =I (r)4 +[{K(A—2)/2}2/4H®] 
« [A/(A — 2) {K(A — 2)/2}2/41P((2A — 4)/A) — 2((A— 1) /A)T'((A—2)/A)]. (35) 


Attention may be drawn to the fact that if the p? factor in (17) is replaced by 
ry ={—a(A—2)/2kT}"T((A—2)/A), (36) 


the term involving it becomes equal to the simpler of the above approximations. 
The value of the ratio r5/p? is 0-67 when AJ is 3, 0-89 when A is 4 and 1-07 when 
A is 6, 

Case 2.—If H<1, U is given by (28) except for velocities well below the 
maximum of the Maxwellian distribution. ‘Taking the limits as before, (27) at 
once yields the approximation 


I(r) = aa (1 —@(H"2)) + H¥?, 
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where ®(H) is again the error integral. It may be noted that the second term 
of (17) can be reduced to this form. A fair approximation (which will be referred 
to as £(r)”) is merely to take 4(r) as unity. 

For intermediate cases numerical methods must be used. Such cases usually 
occur at the large internuclear distances where H and K are identical. To 
facilitate the rapid computation of radiative association coefficients Tables Ia, IB 
and Ic are presented.* They relate respectively to the interaction potentials 
for which A is 3, 4 and 6, and give the values of .4(r) at sufficiently close intervals 
of H (or K)to permit interpolation. The ratios of 4%(r)*, 4(r)*, 4(r)° and 4(r)” 
to the exact function are also displayed. As can be seen from them, the four 
approximations have a remarkably wide range of applicability provided high 
precision is not required. Indeed for many purposes it is sufficient to adopt the 
simple procedure of using 4(r)4 and 4(r)” only, the former inside and the latter 
outside some suitably chosen internuclear distance. 

2.3. If D(r) is the dipole strength of the stabilizing transition, then 


Afr) = “EES! DEW A)— VAN. (38) 


h being Planck’s constant, c the velocity of light and G, the orbital degeneracy 
factor of the final state aust Substitution in (26) gives 


* eae wa |” 


or more conveniently 


I (N)D(X{U yr) — Ur) Pr? dr, (39) 


-o ) 
y= 151 x10 l4gG, | 4(r)D(r){U,(r) — UAr) Pr *} cm.*/sec., (40) 
~ ~O 


where r is now in angstroms, D(r) is in atomic units (ase*), U,(r) and U,(r) are in 
electron volts, and .4(r) (as always) is dimensionless. 

The formula based on a combination of the approximations 4(r)4 and 4(r)? 
can similarly be reduced to 


y= {770 x 10% 37! | Dir) (Ur) — Usr)}(— U(r) dr 


+ 1°34 x 10° "gG, | é D(r){U,(r) — U(r) }r* dr} cm.*/sec., 


r= {—1-1, x 104a/T}*4f(A) sq. angstroms, 
fla) = {(A~2)/2}°T((A—2)/A). (43) 


Here a is some suitable length, such that “(r)“ may be used for r<o and S(r)? 
may be used for r>o; a is the numerical value of the constant appearing in 
(20) when r and U,(r) are in the units specified in the previous paragraph; T is 
in deg. K. The value of f(3) is 1-69, that of f(4) is 1-77, and that of f(6) is 1°71. 

For most transitions D(r) is only known in the neighbourhood of r,, the 
equilibrium nuclear distance. Fortunately, owing to the {U,(r)— U,(r)}* factor, 
the dominant contribution usually arises from this region. Though the constant 


* My thanks are due to Mrs K. Ledsham for assistance in the preparation of these tables 
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TABLE I 


The exact function §(r), and the ratios to it of the approximations $(r)A, 
SI (r)8, F(r)© and J (r)? 
Ia A=3 
H(orK)  S4(r) (AIP (r)  K(r)B (vr) K(v)"/ (rv) F (r)P/ FW) 
0°00 0886 0'00 —2 1°00 1°13 
0°25 1028 0°65 —0'59 1°02 0°97 
0'50 "124 0°67 +or18 1°04 0°89 
o'75 197 0°67 0°43 1°O7 0°84 
1‘00 255 0°67 0°55 1‘1o 0°80 
1°25 303 0°67 0°61 1°13 0°77 
1°50 "342 0°67 0°65 1°16 0°75 
1°75 376 0°67 0°68 1‘1g 0°73 
2°00 "404 0°67 o°’71 1°22 o'71 
3°0 484 0°68 0°76 1°34 0°67 
4°0 "531 0°69 0°79 1°45 0°65 
5°0 563 o'71 0°81 1°56 0°64 
7s ‘610 0°85 1°81 0°62 
10 641 o'75 0°87 2°02 0°61 
15 “689 o°7 0°90 2°37 0°59 
20 728 0°92 65 0°58 
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of proportionality in (20) is rarely well determined, this does not in general cause 
very serious uncertainty. 

2.4. ‘To avoid having to introduce a multiplicity of formulae it has so far been 
assumed that 


(i) —U,(r)>kT for part of the region where the potential has its asymp- 
totic form; 


(ii) U,(r)— U,(r)>kT throughout the region where A(r) is appreciable ; 


(iii) U,(r) <o (except at small r). 


The treatment can readily be generalized to cover other less commonly 
occurring circumstances. If Tis so high that (i) is invalid, then pj, is not given 
by (21); but it may be represented by several such expressions each accurate 
over a limited range. Should (ii) not be satisfied, it is only necessary to take the 
limit €, (in (27)) to equal U,(r) — U,(r) instead of taking it as infinite. The most 
important alternative to (iii) is that U,(r) >o everywhere; and for this case (10), 
with the / term omitted, may be used. Awkward complications however arise 
when U(r) has a positive maximum at large r but is otherwise of the form 
considered earlier; instead of being zero the limit ¢, (in (27)) must be taken as 
equal to the value of U,(r) at the maximum and some suitable algebraical 


expression for pj, must be found. The computations, though not difficult, are 
likely to be tedious. 


3. Applications 
3.1. Because of their possible importance in the study of interstellar clouds, 
calculations were carried out on the rate coefficients y, and y, associated with 
C+ H\(B*>-)— CH(X?M1) + hv (44) 


and 
C* + H\ (AML) CH*(X 15%) + hv (45) 


respectively. _ From the structures of the particles concerned it is apparent, 
cf. (2) and (3), that for the former, 


g=2/18, G,=2, A=6, (46) 
and for the latter, 


g=2/t2,  Ga=1, A=4. (47) 


‘The various potentials can of course be represented by Morse curves modified 
to have the correct asymptotic form.* They were computed from the spectro- 
scopic data given in Herzberg’s recent book (5). The parameter rj (in (42)) was 
rather arbitrarily assigned the value {5 x 10/77} sq. angstroms for CH, and 
{4 x 10% 742} sq. angstroms for CH*+. Accurate information on D(r,){B?X~—X II} 
and D(r,){A 'Il—X 4+} is not available, but for reasons discussed elsewhere (6): 
it seems likely that each is about 3 x 10~* atomic units. Unfortunately, complete 
confidence cannot be placed in this figure. 

On evaluating (41), y, was found to be some 2 x 10718 cm.*/sec. at 100 deg. K..,. 
and to vary as 7” over a considerable range; and y, was found to be also of this 
magnitude, and to be practically temperature independent. 

* It may be remarked that Herzberg and Mundie (4) have found that the A'JI term of the 
molecule BH (which is isoelectronic with the ion CH*) possesses a potential maximum. The 


repulsion implied by such a feature is unlikely to occur when one of the constituent atoms is 


charged; but if it did occur it would of course greatly reduce the radiative association coefficient 
at low temperatures. 
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It must be emphasized that the above results are not necessarily applicable 
under all conditions. One of the factors which makes studies of interstellar 
chemistry unique is that almost every atom is not merely in its lowest term, 
but is in the lowest level of this term. That this is the case for C and C* can readily 
be verified, noting that the Einstein transition probabilities for C(*P,—*P,), 
C(®P,—*Py) and C+(?P3.—?Py.) are 2:8 x1077/sec., 5°83x10-%/sec. and 
2-4 x 10-*/sec. respectively.* The fact that the carbon is mainly C(*P,) and 
C+(? P42) may profoundly influence the radiative association coefficients. Clearly 
the simple weighting ratio g (defined in Section 2.1) is not necessarily correct. 
Account must be taken of the correlations between the molecular and atomic 
levels. ‘There would be no difficulty in doing this if C and C* had strong 7-0/ 
coupling; for the procedure introduced by Mulliken (g) could then be employed. 
This procedure is very simple. It depends essentially on the application of 
Wigner’s rule that the same species of state cannot cross, and it requires only 
knowledge of the order in which the levels of the molecule, and of the separated 
atoms, lie. Provided the #-c coupling is strong it generally gives the correlation 
uniquely. For in the circumstances the c-case of Hund is attained at very moderate 
internuclear separation and hence there is but an unimportant range where, owing 
to the quantum number A being defined, Wigner’s rule is at its least powerful. 
In contrast, if the #-c% coupling is weak the range for which is a good quantum 
number is extremely wide, and as a consequence Wigner’s rule (which allows 
states of different A to cross) is insufficient to determine the correlation. ‘The 
coupling in both C and C* is very feeble indeed, as evidenced by the fact that 
the extreme separation of the ground term levels is only 43 cm.~! for the atom 
and 64 cm.! for the ion. Mulliken’s method cannot therefore be used. 
A rigorous alternative method is not easy to develop, but the following is 
suggested as plausible. 

Consider the effect of the axial electric field on a pair of atoms as they approach 
each other. If the #-c% coupling is sufficiently weak it is broken at large inter- 
nuclear distances where the system can still be described by assigning quantum 
numbers to the individual atoms. ‘These quantum numbers can be found from 
the Stark pattern (10); and knowing them it is possible to place restrictions on 
the correlations between the atomic and molecular levels. Thus in an electric 
field C(?P,) probably becomes C(*X), and C*(?P,).) probably becomes C’(*Z), 
so that addition of H(#S) to the distorted atom can presumably only yield CH(?2) 
or CH(*X), and its addition to the distorted ion can presumably only yield CH *(#Z) 
or CH*(8X).¢ Now in the formation of CH and CH* by radiative association 
the initial terms involved are B*X~ and A'Il respectively. Provided the move- 
ment is adiabatic it would Seem that the chance that C(*P,) and H(?S) come 
together along the former is }, and the chance that C‘(?P,,.) and H(?S) come 
together along the latter is zero. Using these values for the factor g appearing 
in the original formula (instead of the simple weighting ratios adopted previously) 
gives finally y3=6 x 10-8 cm. /sec., (48) 


Ya=0. (49) 


* The first two values are taken from the tables of Pasternack (7); the third value was calculated 
from the formula given by Shortley (8). 

+ It is satisfactory to note that similar arguments applied to CdH lead to results in agreement 
with the observation of Svensson (11) that B*Z is correlated with Cd(*P,)—H(?S) and that A*/7, , 
is correlated with Cd(*P,)—H(S). As Mulliken (9) has pointed out, this correlation is not given 
by the other scheme. 
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It may be remarked that if the #-c/ coupling were strong the correlation 
would be such that both coefficients would be zero. As the two extremes thus 
agree in the case of CH* it seems likely that the true coefficient is indeed very 
small (unless the motion is not closely adiabatic). The position regarding CH 
is perhaps more doubtful. The coupling is certainly extremely weak, and one 
would therefore, for the reasons given, expect the coefficient to be comparatively 
large. Nevertheless it must be borne in mind that the branch of molecular 
theory concerned is still almost completely unexplored, and consequently the 
correlation scheme proposed should be treated with some reserve. In particular 
it is desirable to investigate the probability of non-adiabatic collisions*, and to 
ascertain the reason for the predicted crossing of the CH, X 7II1,;. and B*=- 
potential energy curves. Mulliken (12) suggests that the final analysis may 
require detailed calculations, taking into account the interactions between the 
rotational, vibrational and electronic motions. 

Bates and Spitzer (6) have studied the equilibrium in interstellar clouds, using 
both the sets of coefficients. ‘They find that neither set is sufficiently large to 
explain the observational data if the conditions in the clouds are as generally 
believed (13). ‘Though correlation schemes could be visualized which would 
lead to a decrease in the coefficients, none could lead to a significant increase. 
The discreparcy would thus seem real. 

3-2. Knowledge of the rate coefficient, y,, associated with 


N++4N|(B*E)—> Ni(X2D3) + hv (50) 


is required in certain problems arising in the study of the Earth’s upper atmosphere. 
In the case of this process 


g=2/36, G,=1, A=4 


(cf. (2) and (3)). As usual r5 (of (42)) is poorly determined. It was taken to be 
{5 x 10?/71"} sq. angstroms. For the applications envisaged the conditions are 
such that consideration need not be given to the detailed correlations. There 
is still controversy regarding the dissociation energy of Nj. To minimize the 
possible error the mean of the values advocated by Gaydon (3) and by Herzberg (4) 
was adopted in constructing the relevant Morse curves. According to the 
calculation of Shull (14) D(r,){B*2X —X?Z7} is about 1-5 atomic units. 

Taking the various quantities appearing in (41) to be as indicated above, it 
was found that over a wide temperature range (up to at least 2000 deg. K.) 


Y3==3 x 107!’ cm.3/sec. (52) 


The suggestion has been made (15) that process (50) is responsible for the 
presence of Nj ions in high-altitude sunlit aurorae. It is apparent however that 
the magnitude of y, is much too small for this to be accepted. For even if the 
concentration of N atoms near the 1000 km. level were as high as 108/cm.* the 
lifetime of an N* ion towards radiative association would be of the order of 
10 years. 

The extent to which the oxygen in the upper atmosphere is dissociated depends 
partially on the important reaction 


O+0O)|(A®Zf)—> OX ®Z>) + hv. (53) 





* If this probability is actually large (as might be the case since the energy differences are 
minute) the normal weighting ratio g should be used. 
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Unfortunately the rate at which this proceeds cannot be determined with any 
precision at present as the dipole moment of the transition is not known reliably. 
It is interesting to observe that the photon yield P may be high. If y, is the 
coefficient involved, n(O) is the concentration of atomic oxygen in the region 
where it is most abundant, and H is the local scale height, then 


<n(O)? Hy,/cm.? column sec. (54) 


Since n(O) is about 10!/cm.’, and H is about 10°cm. (16), P is thus of order 
10**y,/cm. column sec., and must therefore be very considerable unless y, is 
extremely small.* It would be of great value to measure the intensity of this 
continuum. 

3.3. Collision processes involving hydrogen only are of special importance in 
astrophysics and consequently it was considered worth evaluating the rate 
coefficient, y;, associated with 


H+ H*\ (2X) Hy(X*X7) + hv. (55) 


Additional interest in this reaction is provided by the fact that it is an example 
of the class in which the approach is along a repulsive potential energy curve. 
A welcome feature of the calculations is that they can be carried out with some 
precision (for essentially the only approximation is the neglect of the wave nature 
of the nuclei). 


TABLE II 
Radiative collisions between hydrogen atoms and tons 


‘Temperature Rate coefficient for Sum of free-bound and = Mean energy 
(T) Hy formation free-free rate coefficients of emitted 
(ys) (ys+¥5) photons 
(deg. K.) (cm.3/sec.) (cm.3/sec.) (e.V.) 
500 :3xX26"" ‘ so" 23 
1000 $°2X 10°" K167>° “41 
2000 9X 107" tac’ "72 
4000 6-2 x 10717 x 19~1¢ 
8000 17x 107° 19~'6 
16000 3°53 x t0"™ K 107 
32000 5°6 x 10716 K 10-% 
64000 6:0 x 190716 K 10-15 


CwWwrena 
-$nonorss 
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In a recent quantal study of the two-centre problem Bates, Ledsham, 
and Stewart (17) have derived the potentials concerned, and Bates (18) has used 
the wave functions made available from this work to compute the Einstein 
coefficient at various internuclear distances. All the quantities in (11) are thus 
known accurately. The necessary integration was performed by numerical 
methods and hence from (10) (in which the weighting ratio g is here }) y; was 
found. ‘Table II shows the results. It will be noted that the coefficient increases 
rapidly with the temperature, as would be expected. 

All radiative collisions do not of course lead to the formation of Hj —after 
some the H and H* particles can move apart again (with decreased kinetic energy). 

* Thus P would be some 1o'?/cm.? column/sec. if, for example, y, were 10~*® cm.*/sec. As a 


comparison it may be noted that the number of photons emitted in the strongest line of the 
normal airglow is only about 10°/cm.? column/sec. 
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The sum of the free-bound coefficients y, and the free-free coefficient y; can be 
obtained by putting A(x.) in (10) equal to unity. For completeness this sum 
is also displayed in Table II together with the mean photon energy emitted per 
radiative collision.* 
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Summary 


Photoelectric observations of the system { Phoe are described. The system, 
which has not hitherto been recognized as light variable, undergoes total 
eclipses. Uniform and partially darkened solutions are given. The orbital 
eccentricity is determined solely from the photometric results. - It is notably 
less than the figure derived from the spectroscopic data. Certain other 
discrepancies in the solution are discussed. These can be resolved by 
assigning to the components temperatures which are much lower than would 
be expected from the spectral types. Absolute dimensions are estimated. 





1. Introduction.—Recently some photoelectric measures were made of 
certain southern stars which were known to be spectroscopic binary systems but 
which, as far as could be ascertained, had not been recognized as light variables. 
Amongst these was the system ¢ Phoenicis (« =o1" 04™-2, 8 = —55° 47’ (1900), 
B8 + B8, 4™-13) which had been reported as a spectroscopic binary by R. E. 
Wilson from Lick observations at Santiago. Later A. Colacevich* assigned a 
period of 19-66958 to the system and deduced a set of elements. ‘Iwo spectra 


were reported, but the lines of the brighter were broad and those of the fainter 
were found on but five plates and then only with difficulty. Because of uncer- 
tainty in the measurements Colacevich regarded the elements as provisional. 
‘The photoelectric observations described here have shown the star to be a well- 
marked eclipsing system, one of the few to be discovered photoelectrically. 
The system ¢ Phoe is a visual triple. SDS Rmkz lists a pair 4™-2 + 8™ separated 
by 6” and effectively stationary over the last century. Rossitert+ observed a 
third star 7™-2 separated by 0"-6 from the 4™-2 component (Rst No. 1205). 
The present observations refer to the brighter component and any perturbations 
introduced by the fainter companions have been disregarded. 

2. Photometric technique.—The star was observed with the aid of a photo- 
multiplier tube attached to the 6-inch Farnham refractor. ‘The photometric 
methods were similar to those used in the observation of S$ Antliae{ in that the 
output meter of the amplifier was used as an indicating instrument only, the 
actual measurement being made by a laboratory potentiometer. A sketch of 
the circuit is shown in Fig 1. The 931A multiplier tube was modified so that 
the output, instead of being through a pin close to dynode pins at a high voltage, 
was led out through an earthed guard-ring in the side of the socket. ‘The high- 
voltage supply was obtained from a controlled rectifier set. ‘This offered a 
convenient method of adjusting the sensitivity of the photometer. ‘The amplifier 
tube was run much below its rated voltages so that the grid current was negligible 
in comparison with the photo-current. 


* A. Colacevich, P.A.S.P., 47, 84, 1935. 
+ R. A. Rossiter, Mem. R.A.S., 65, 40, 1933- 
} A. R. Hogg and P. W. A. Bowe, M.N., 110, 373, 1950. 









































Fic. 1.—Diagram of photo-multiplier circuit. 
V , =600-g00 v. Vy =45 Vv. V3=4'4V. Vi=1'5v. 
V;=22¢ v. R, =50 megohms Rz= 10,000 ohms. 
C,=0'1 pfd. MP=measuring potentiometer. 


TABLE I 
Normal Magnitude Differences « Scul-€ Phoe 


The normals are the means of ten observations arranged in increasing order of phase 
for a period of 14-6698 with zero phase adopted at J.D. 2432499000 


Phase Am Phase Am Phase Am Phase Am 
d m d m d m d m 
0°036 0°303 0°376 0°304 0°976 O°'112 mss 0°279 
“076 +302 “416 *307 -985 ‘O51 "165 "204 
"102 ‘280 *434 ‘296 ‘000 —0'029 “176 "294 
"119 "226 "454 *308 ‘009 —0'094 "194 *300 
"135 “198 *489 "312 “o18 —O'I§4 "212 *297 
“146 "153 "510 "302 "026 —0*167 °232 "299 
“159 ‘217 °527 *313 "030 177 ‘278 "294 
"174 ‘077 "547 "312 "O41 188 "317 "299 
"188 "057 "574 "310 ‘O51 170 "369 278 
"203 056 ‘620 *318 ‘060 175 "414 295 
"215 "044 681 *300 ‘O71 139 "450 314 
*230 064 -766 *320 ‘o88 "024 “485 “yy 
"244 "155 “S18 "299 *100 "O55 "530 *303 
*256 ‘IS! 853 303 "109 116 "550 "330 
°275 "215 “867 *298 ‘119 *167 "567 "310 
"290 "256 ‘882 "304 *126 “198 *579 *310 
*305 282 *g00 *302 133 +228 591 302 
*323 +296 "929 *300 “140 *263 606 *308 
0°348 0°295 0960 = 0203 1°47 264 628 "295 
1648 o-291* 


000000000 


The comparison star used was « Sculptoris (« =o" 53™-8, 5 = —29°54’ (1900), 
4™-39 B5), which although rather distant from the variable was of comparable 
brightness and spectral type. A total of 772 comparisons on 63 nights between 
1947 November and 1949 February was obtained. A comparison usually 
consisted of one setting on the variable, interpolated between two settings on 


* Four observations only. 
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the comparison star. Five readings of the potentiometer were made for each 
setting. Because of the small size and convenience of the instrument a complete 
setting (finding and reading) could be done in three minutes or less. Corrections 
for atmospheric extinction were determined and applied as described elsewhere. * 
The observations were assembled on a period of 14-6698, and the normals shown 
in Table I were derived. Each normal represents the mean of ten observations 
as arranged in increasing order of phase. 

3. Period.—The period adopted in Table I (14-6698) differs slightly from the 
value used by Colacevich (14-66958). The adopted value was derived from the 
earlier portion of the photometric results and in its turn differed slightly from an 
improved value obtained from all the present observations. The following 
heliocentric times of minima were observed, viz. :— 

Primary Minima J.D. 2432772217 
797°263 
818-978 
849°025 
Secondary Minima 813-128 
889-949 
899-963 
These results gave 
J.D. Epoch of principal minimum = 2432667-012 4- 1-66990 E. 

4. Rectified light curve.—A plot of the normals (not shown) indicated the 
presence of elliptic components. Accordingly the points outside of the eclipse 
were formed into ‘“‘super-normals’’ of 20 observations each and fitted by least 
squares to the usual expression 


L = Ay+ A,cos@ + A, cos 20. 


(ZL =light intensity, A,=reflection coefficient, A, ellipticity coefficient.) The 
results obtained were as follows 


Ay = +0°9861 +0-0009 (p.e.), 
A,= +0°0005+0°0014 (p.e.), 
A, = —0:0107 +0-0016  (p.e.). 
Thus the reflection effect (A,) was effectively absent and the ellipticity (A,) was 
small but statistically significant. The observed points were rectified for ellipticity 
by division by (A)+ A, cos 20) in the usual manner. ‘The rectified observations. 
were then folded (a) about the primary minimum and (4) about the secondary 
minimum. ‘The folded poirts within the eclipses were each derived from two 
sets of sub-normals each of five observations taken at approximately equal times 
before and after the minimum. The points outside the eclipses were derived 
in a similar manner from pairs of normals given in Table I. ‘Thus for the folded 
rectified curves the final points within the eclipses each represent 10 observations 
(with one exception at secondary minimum) and those outside the eclipses each 
represent 20 observations.. These folded rectified curves (Fig. 2) strongly 
suggested a total eclipse. The flat bottom to the secondary curve indicated an 
occultation eclipse, whilst the slight curvature of the bottom of the primary 
minimum indicated a transit eclipse (perhaps partial or perhaps total) across a 
darkened disk. 
* A. R. Hogg, M.N., 106, 294, 1946. 
22* 
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5. Uniform solution—A solution for uniformly illuminated disks was 
obtained by the methods due to Russell and Shapley.* The rectified depths of 
the minima from a curve drawn for total eclipses were primary =0°3525, 
secondary = 0-1955 (in light-intensity units). The Russell-Shapley equation (j), 
with the hypothesis that the principal minimum was an occultation, then gave, 
for a total eclipse, k (the ratio of the radii)=0-550. The hypothesis that the 
principal minimum was a transit gave k=0-662 for a total eclipse. The light 
curve for primary minimum, when analysed with the aid of Russell’s -functions 
and the assumption of a total eclipse, gave k=0-647. Rather more weight was 
given to this latter value as the individual values of k from the ¥-functiong agreed 
quite well, and the value of k=0-65 was adopted. 

The light curve showed unequal spacing of the minima, thus indicating the 
existence of an eccentric orbit. The average time of occurrence of each minimum 
was derived from a curve drawn through the normals of both minima. Eight 
equidistant light levels were marked on the curve by a series of parallel lines, 
each one cutting the ascending and descending branches of a minimum. The 
mid-point of each parallel between these points of intersection gave a value for 
the time of minimum based on a symmetrical curve. The means of these points 


gave Phase of primary minimum =o0"-000 +0%004 =, 


mi secondary __,, =0"-8257 + 0-009 =f. 
As h=ecosw=(t,—t,— P/2)/ P(1 + cosec*i) = —0°0027 + 0-0003 (p.e.), some 
eccentricity was indicated. (Cosec?i was taken as 1-0, an assumption justified 
by later results.) This result of a finite value for A indicated that the value of w 
would ditfer from the spectroscopically determined value of w=go°. Accordingly 
an attempt was made to derive the value of g =e sinw from the light curve by the 
separate analysis of the primary and secondary minima using a common value of 
k as described by Russell and Shapley. Kopal{ has shown that if sin®@ is 
plotted against [1 + kp(k, «)]? a straight line will result, from the slope and intercept 
of which values of r, and i may be obtained. ‘This method is applicable directly 
to normals. It has been applied to the present observations with the necessary 
modification to take account of the ellipticity of the components, i. e. the division 
of sin? @ by 1—2cos*@, where the photometric ellipticity constant was taken as 
z=4A,. Fig. 3 shows the graphs for primary and secondary minima. The 
points graphed are the actual normals, not smoothed values. The lines shown 
are the regression lines of the term in p on the term in 6. ‘The results were as 
follows : 
Primary Minimum Secondary Minimum 

Semi-major axis, larger star a’ =0°250 a” =0'261 

Orbital inclination 7 coti’ =0-0508 coti” =0:0566 

7» = (2 + cos Z)§ (gn)' =0-0218 (gn)” =0-0284 
Using Russell’s equations, aj = a;(1 + 2gn) and coti” =coti'(1 + 4gy), it was found 
that the mean value of gy (see above) =0-025, from which g=esinw=0-026. 
Combination of this with h gave, for the uniform solution, 


e =0°026 w =96°-0. 





*H. N. Russell and H. Shapley, Ap. 7., 36, 404, 1912. 

t Loc. cit., p. 54. 

t Z. Kopal, Ap. 7., 94, 149, 1941. 

§ ¢= (6 - 6”), where & and 6” are values of @ at beginning or end of the eclipse. 
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For small values of gy the true values of a, and coti approach the means of the 
values derived from the minima separately, which are shown in Table II, along 
with other quantities derived by the usual methods. 
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Fic. 3.—Characteristic diagrams for primary and secondary minima. Uniform solution. 


Kopal’s method. 


Taste II 
Solution for Ellipsoidal Stars in Eccentric Orbit 
Period = 14-66990 
Darkening Coefficient 

Semi-major axis, relative orbit 
Semi-major axis (equatorial), larger star 
Semi-minor axis (equatorial), larger star 
Semi-axis (polar), larger star 
Semi-major axis (equatorial), smaller star 
Semi-minor axis (equatorial), smaller star 
Semi-axis (polar), smaller star 
Axial ratio (a,/a,; = b,/b, = ¢2/c,) 
Light of larger star 
Light of smaller star , O°195 
Equatorial ellipticity coefficient = e? sin? 7 0°040 
Inclination of orbit (corrected for polar flattening) 7 87°-o 
Eccentricity of orbit ? 0-026 0'028 
Longitude of periastron 96°-o 95 °6 
Fractional light loss Xo 1‘0o 1'004 
Ratio of surface brightness J,JJ, 1°80 1°84 


6. Darkened solutions.—Trials were made to obtain completely darkened 
circular solutions for the two minima separately on the basis of the secondary 
minimum being an occultation. It was found that satisfactory solutions with 
a common value of k as required by the eccentric orbit theory could not be 
obtained. Using Merrill’s tables* other values of the darkening coefficient x 

* J. E. Merrill, Princeton Contr., No. 23, 52, 1950. I am much indebted to Dr J. E. Merrill 


for promptly despatching a copy of these most valuable tables, which greatly helped in the darkened 
solution. 
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were tried. The most suitable was x =0-2 and although this was rather smaller 
than the value normally adopted for the spectral type of the components it was 
used for further work. The value of the ellipticity constant z appropriate to 
this value of x was derived from the constant of the (cos 20) term in the analysis 
of the inter-eclipse points, using the expression due to Russell*, 


z= ¥ (4,003) 


where the value «/C, representing the ratio of the geometric to the photometric 
ellipticities, was taken as 2-7 from Russell’s data. This method gave *z =0-038 
which was used for the phase rectification in the usual manner.f Assuming a 
total occultation at secondary minimum and using Merrill’s tables of * in 
conjunction with Russell’s equations{ gave a value of k=0-625 from the light 
curve on the basis of rectified circular elements. Circular elements were next 
derived for the primary eclipse with the same value of k =0-625, using Merrill’s 
tables of 2" and eight points on the light curve. This gave ap” =1-004 with a 
value of p= —1-057. The elements, determined in the same manner as for 
the uniform solution, together with the values of g, were as follows :— 
Primary Minimum Secondary Minimum 

Semi-major axis of larger star a’ =0°262 a” =0°272 

Inclination 7 coti’ =0-0878 coti” =0*1005 

&7 (gn)' =0-019 (gn)” =0-028 
Again using a mean value of gy and combining the value of g with the value of h 
from the spacing of the minima gave e=0-028 and w=95°-6. As in the uniform 
solution, the actual values of a, and coti for this darkened eccentric solution were 
obtained as the mean of the circular elements for the two separate minima, 
viz. €y=0°267 andi=84°-7. (Table IT.) 

A characteristic diagram, using the observed normals for the partially darkened 
solution, similar to that for the uniform solution, is shown in Fig. 4. For the 
partially darkened solution the more convenient functions from Merrill’s tables 
have been used instead of the p-functions. At the time the uniform solution 
was made Merrill’s tables were not available. Fig. 4 also displays the differences 
between primary and secondary minima indicative of an eccentric orbit. 

Both the uniform and the partially darkened solutions gave satisfactory fits 
to the primary and secondary minima, but the bottom of the primary minimum 
was better represented by the partially darkened solution. The line in Fig. 2 
represents the solution for the darkening coefficient x=0-2. For the paired 
sub-normals in the primary eclipse the s.d. of a single point was 0-003 light units ; 
for the secondary eclipse the s.d. was 0-004 light units. The “super normals” 
between eclipses gave a s.d.=0-005 units. No deviation exceeded 0-010 
light units. 

7. Ellipticity of components.—The foregoing solutions are based on the 
components being similar and similarly situated ellipsoids with their longest 
axes coinciding with the line joining their centres. This physically unlikely 
model may be brought nearer reality by assuming a degree of polar flattening. 
The polar flattening cannot be derived from the photometric results but has been 
introduced into the solution by Russell{ who gives <? =z cosec”i, where ¢ is the 

*H. N. Russell, Ap. 7., 102, 1, 1945. 


+t H. N. Russell and H. Shapley, Ap. 7., 36, 407, 1912, equations c’ and f’ 
$H. N. Russell, Ap. 7., 36, 62, 1912. 
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eccentricity of the meridian section such that b = a(1 —e*)"* and from the dynamical 
theory for small amounts of «,c =b—}(a—6). Here a, b and c are the axes of the 
ellipsoid in descending order of magnitude. ‘The assumption of polar flattening 
requires that the value of i should be altered so that cot i;=(c/b) cott,, where the 
subscripts e and f refer to the elliptical and the flattened conditions respectively. 
Values derived on this basis for both uniform and partly darkened solutions are 


shown in Table II. 
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Fic. 4.—-Characteristic diagrams for primary and secondary mimima. 
Partially darkened solution (x =0-2). 


Certain relations should exist between the term A, in the inter-eclipse light 
curve and the elements. Results based on dynamical theory, combined with 
a gravity- and limb-darkening expression, led Russell* to the conclusion that the 
ellipticity coefficient should be given by 


A, = [ e073 (7 r+ ee r) -K(C,L,- Cals) | sin? i, 


where, in addition to the usual symbols, E, and E, are the luminous efficiencies 
at the wave-length measured for large and small stars respectively, 7, and rz 
the mean radii—(abc)"* and C, and C, the computed photometric ellipticities. 
Using the derived values of r, L andi, together with results tabulated by Russell, 
and taking E, = E, as the components are of the same spectral type, the computed 
value of A, is found to be = —6-0099. This is in agreement with the observed 
value of —0-0107 + 0-0016. Thus the computed values of L, r and i are consistent 
with the observed values of the ellipticity coefficient. 

8. Absolute dimensions, luminosity, density, etc.—The combination of the 
photometric results with the spectroscopic figures of Colacevich (l.c.) gave the 
quantities shown in Table III, which refers to the x=0-2 solution. The stars 
appear to be main-sequence types. Colacevich, using the mass-luminosity 
relation, derived absolute visual magnitudes of 1-6 and 3-0 for the larger and 

* H. N. Russell, Ap. F., 104, 153, 1946. 
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smaller components respectively. ‘This difference of 1™-4 between the com- 
ponents is in good agreement with the present work, which gives (from the ratio- 
L,/L,) a difference of 1™-5. Other comparisons with the spectroscopic results. 
disclose certain difficulties. ‘These are as follows :— 
(1) The value of J,/J, (and T,,/7,,) differs from that expected. 
(2) Computed masses differ from the spectroscopic figures. 
(3) The density of the faint component is several times greater than would 
be expected from its spectral type. 
(4) The absolute magnitude differs from that to be expected from spectro- 
scopic parallax results. 
TABLE III 
Absolute Values from Spectroscopic Orbit, Darkened Solution 
and Spectroscopic Parallax (0”:014) 
Star 1 Star 2 
Semi-major axis, relative orbit km.» 10% 2°47 4°07 
Semi-major axis of star km. 108 = =61°74 =~ 1°09 
Mean radius, (abc)! Solar units 2°5 1°5 
Mass Solar units 2°8 2°0 
Density Solar units 0-18 0°54 
Absolute magnitude o'2 1°79 
‘Lhe ratio of the surface brightness of the components, J,/J,, was derived from 
a comparison of the depths of the rectified minima as 1-8; but as the two com- 
ponents have the same spectral type a value near unity might be expected. This 
discrepancy was also demonstrated by making use of an expression for computing 
the photographic absolute magnitude from the radius and effective temperature 
of a black-body star.* The effective wave-length of the photographic system 
is taken by these authors as 4250 A. and is probably near the effective wave-length 
of the photomultiplier system. From the above-mentioned expression it is 
easy to show that , ee ° kL, 
ma — a =0°52 xX I0-° lo ; 
ae ema: © 
Adopting a value of 7',,=12,500 deg. as appropriate to a BS star gave 
T,.=10,700 deg. ‘The effective temperature of the smaller star would thus be 
equal to that assigned to an AI type. 

For computing masses, Gaposchkin*t gives an equation combining Kepler’s 
third law, the Stefan-Boltzmann expression and the mass-luminosity relation. 
Using this equation it was found, if values of 7, = 12,500 deg. and m,/m,=0-71 
(spectroscopic value) were adopted, that m,=5:3 and m,=3°8 in solar units. 
These values are about double those obtained spectroscopically. 

The third discrepancy between the photometric and spectroscopic observations. 
lies in the derived densities. ‘Thus p, determined from the spectroscopic value 
of the masses and orbital dimensions combined with the photometric values for 
the radii and 7, is 0-18 solar units, i.e. about equal to the density to be expected 
in a main-sequence star of the spectral type BS. On the other hand, a similar 
estimate of p, gave 0°54 which might be attributed to a spectral type A2 or later. 
This does not agree with the spectral observation of the two components being 
both of the same type, viz. B8. Also the values of the density derived from 
Russell’s{ equation 52, which depends on Kepler’s third law and involves. 
luminosity only as a correcting term, were p,=0-18 and p, =0°43. 

*H. N. Russell, R. S. Dugan and J. Q. Stewart, Astronomy, I, p. 732, 1927. 


+ S. Gaposchkin, Proc. Amer. Phil. Soc., 79, 329, 1938. 
tH. N. Russell, Ap. 7., 36, 74, 1912. 
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The fourth discrepancy arises when the absolute magnitude, as determined 
by the classification of objective prism spectrograms, is used with Kuiper’s* 
bolometric corrections and empirical mass-luminosity relation and the present 
photometric value of L,/L,, to obtain separate masses for the components. 
Unpublished results by Miss Woods resulting from a re-classification of spectro- 
grams taken by Rimmerf were used. These results gave an absolute magnitude 
of 0-0 (pv) for the complete system, which with the photometric ratio L,/L, 
gave figures of 0-2 and 1-7 for the large and small components respectively. 
These, with Kuiper’s results, gave m,=4:1 and m,=2°8, in not over-good 
agreement with the results from the orbit. ‘These discrepancies are all associated 
with the various assumptions made about the effective temperature. If tempera- 
tures of 8250 and 7450 deg. could be assumed then, approximately, 

(1) The equation derived from Russell, Dugan and Stewart would be 
satisfied. 

(2) Gaposchkin’s equation would give masses in agreement with the 
spectroscopically observed figures. 

(3) The observed densities would be more in keeping with those to be 
expected from the spectral classes associated with the assumed tempera- 
tures. 

(4) A different bolometric correction is applicable and accordingly the 
values of the mass from the objective prism results become m, = 3:2 and 
m, = 2°1 solar units in quite satisfactory agreement with the orbital results. 


Such an assumption however would be quite out of keeping with the assigned 
spectral types and their corresponding temperatures. From Colacevich’s 
remarks there might be some doubt about the spectral type of the fainter com- 
ponent. Granting that the larger star was of the type given by Colacevich (B8) 
and assuming that the fainter star was about half a class later would overcome 
the J,/J, difficulty but would not account for the masses and densities. { 

Conclusion.—The system of { Phoe is of considerable interest and worthy of 
further study, not only toimprove the existing results but also to ascertain whether, 
with the passage of time, any perturbations due to the visual components might be 
detected. Both the direct spectroscopic parallax (0’-014) and the parallax 
derived from the orbit masses with the usual temperatures (0”-033) are large 
enough to suggest that a fair trigonometric parallax might be obtained. 
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* G. Kuiper, Ap. 7., 88, 446 and 489, 1938. 

+ W. B. Rimmer, Mem. Comm. Obs., 2, 1, 1930. 

} Subsequent to writing this section the author has had a conversation with Dr Kron who 
states that he and Dr Herbig are of the opinion that the Lick plates show the brighter spectrum to 


be B7 and that the fainter spectrum, whilst not of high quality, is definitely later. The ratio 
L/L, could be 4. 
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Summary 


The magnitudes of bright stars have been determined, using the band 
4500-4600A. Dispersion was effected with an objective prism spectrograph 
and the light intensities were measured with an electron multiplier in a null 
circuit. 

Seven of the stars, forming a ring around the sky between 6= — 20° and 

70°, were treated as primary stars and each one was compared with the 
six other stars in the group. The standard deviation of a magnitude in this 
group was estimated as +0™-o02 (internal error of the system). The 
magnitudes for other stars were determined by comparison with stars of the 
primary group, using either two stars from this group (secondary) or a single 
star (tertiary). Comparisons have been made with other scales. 





1. Introduction.—In defining a scale of stellar brightness there is a choice of 
a range of alternatives. At one extreme the scale may aim to express the ratio, 
from star to star, of the emission in the continuous spectrum at a particular 
wave-length. At the other extreme it may be based on the response of a particular 
receptor (e.g. an eye) to a relatively wide spectral region with consequent 
peculiarities due to well-known instrumental, atmospheric and physiological 
effects. ‘There are practical advantages in using monochromatic light as the 
basis of any system. Thus (a) the effective wave-length of any measurement is 
independent of the colours of the sources involved, (6) greater precision is 
introduced into the necessary determination of atmospheric absorption coefficients 
which are then independent of colour of the source and zenith distance of the 
measured object, (c) any corrections for drifts in the colour response of the 
apparatus, e.g. change in colour of mirror, are unnecessary, (d) there is greater 
certainty in converting accurately to another magnitude system of defined wave- 
length, and (e) with a monochromatic system it is possible to specify, quite 
rationally, an absolute energy flux value (e.g. a specified number of watts/cm.”*) 
as the zero of the magnitude scale. There are however practical disadvantages 
also, including the risk, as the band width is greatly decreased, of encountering 
appreciable absorption (or emission) lines varying from spectral type to spectral 
type. This could render the results useless, from a continuous spectrum view- 
point, without the prohibitive labour of measuring the absorption in large numbers 
of stars. 

As a compromise between these factors the scale may be based on a narrow 
but not infinitesimal band of wave-length located in a region of the spectrum 
relatively free from heavy stellar absorption lines. ‘This compromise has been 
effected in the present work where the chosen band, between A 4500 and A 4600 A. 
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is free from any very pronounced lines, does not have an excessively great atmo- 
spheric extinction and is in a region covered by the more sensitive of the 
present-day photoelectric equipment. 

The results from such a process could differ from theoretical monochromatic 
magnitudes of the continuous spectrum based on some black-body concept in two 
ways, Viz. 


(i) purely on account of finite band width, 
(ii) on account of stellar absorption lines. 


Effect (i) was investigated theoretically for black bodies. It was found that the 
magnitude difference between two stars of spectral types A and K would vary by 
less than 0™-001 when measured with band widths of 1 A. and 100A. centred at 
4550A. Effect (ii) has recently been discussed by N. Milford (1), who estimates 
that it may be very considerable, especially in later spectral types. Thus in 
comparing « Boot (K2 pec) with B Orio (B8) the correction would exceed 0™-3 
for the band used in the present work. 

From a practical point of view the use of a very narrow band width (say 1 A.), 
especially with stars having many absorption lines, could give results which would 
be quite sensitive to minor variations of the band width used. The use of a 
medium band width (say 100A.) could be expected to give more reproducible 
results, because of an averaging effect between any absorption lines at the longer 
and shorter wave-length limits of the band. It is of course highly desirable that 
the cut-offs of the band should be square, otherwise the effective wave-length 
could alter with the colour of the source. 

The figures presented here refer to the emission (uncorrected for stellar 
absorption) averaged over a much narrower band than is commonly employed in 
magnitude measurements. ‘They may, particularly in the case of later types, 
deviate very considerably from monochromatic magnitudes referred to some 
theoretical continuous spectrum. The bulk of the present work deals however 
with B types where such correction appears to be small. 

2. Equipment.—The present measures were made by an electron multiplier 
photo-tube attached to an objective prism spectrograph on the Reynolds reflector 
(30-inch). The objective prism spectrograph has been described by Gascoigne (2) 
and is shown in Fig. 1. A defining slit in the focal plane of this instrument was 
so placed as to pass radiation between A 4500 and A 4600 when the line HB was 
centred on a locating cross-wire mounted in an adjoining slit. So as to render 
HB more easily visible the locating wire was viewed through a special eyepiece 
containing a cylindrical lens. The light passing through the defining slit was 
collected by a lens (Fabry system) and fell on to the cathode of the photo-tube. 

In setting up the apparatus particular attention was paid to ensuring that the 
expected amount of light entered the equipment. Thus with photo-cell removed, 
the eye was placed near the focal plane of the spectrograph when directed to a 
star and images of the telescope mirrors were examined to determine that all were 
centred and full of light. Further, to minimize the possible effect of any flexure in 
the telescope a circular diaphragm of diameter about } inch less than the aperture 
of the collimated beam was placed in front of the spectrograph prism. Slight 
flexure was present but tests indicated that its effects were negligible. The reflector 
and spectrograph having been aligned on a star showing HB prominently, the 
locating slit was adjusted by a screw until HB appeared in the eyepiece as being 
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on the cross-wire. ‘The defining slit then passed the required band of radiation. 
The guiding telescope (Fig. 1) was then moved bodily until the star image was seen 
on its cross-wires. ‘These cross-wires then provided a fiducial mark when 
measuring stars in which Hf could not be discerned. 


| 


Diagonal 


Guiding Telescope 























Fic. 1.—The optical train. 


For certain comparisons a diaphragm or a filter was employed. ‘The dia- 
phragm of the Maltese-cross type was inserted at CC (Fig. 1). ‘The two filters 
used were of glass coated by evaporation with a /Vi-Cr alloy (3). The filters 
were spectrally neutral to a high degree but did not prove sufficiently stable for 
our purpose despite having been baked to 200 deg.C. in manufacture. One 
filter (2™) was used in the parallel beam at position CC; the other (05) was 
placed in the photometer mounting. Filter transmissions were determined at 
the telescope, using stars near the zenith on good nights. ‘The Maltese-cross 
diaphragm was not suitable for the greatest precision because of possible local 
changes in the reflectivity of various.parts of the mirror surface. Also alterations 
in its transmission could be expected with changes in the spectrograph setting 
when the instrument was replaced on the telescope after making way for other 
programmes. The evaporated film filters were better in this respect, but 
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because of their instability, the later stages of the work were arranged so that a 
knowledge of filter transmission was not required, i.e. full reliance was placed 
on the response of the cell being linear. Filters or diaphragms were then 
used only to bring the light intensity of the brighter star of a comparison to a 
convenient value and were retained in position for the measurement of the other 
star. 

The photometer proper consisted of a selected 931A electron multiplier. 
The anode lead of the multiplier was detached from its pin and brought out at 
the side of the base through an amber insulator surrounded by an earthed guard- 
ring to minimize leakage effects. The output from the multiplier was passed 
through a 10 megohm resistor and the voltage drop over this resistor was measured 
by a single-valve amplifier using a null method. The circuit is shown in Fig. 2. 
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Fic. 2.—Photometer circuit. 
The multiplier tube 931A was operated at 90 v./stage. MP=measuring potentiometer. 


When no current flowed through the photo-cell, the grid of the valve was at or 
near earth potential and the anode current was adjusted by altering the potential 
of the cathode with respect to earth until the microammeter read 25 microamps. 
(The valve was greatly under-run to ensure a grid current below 10~! amp.) 
When the multiplier was illuminated the anode current of the valve was thereby 
decreased. It was then restored to 25 microamps. by altering the potential of 
the cathode, using a potentiometer calibrated in millivolts. Thus the light 
intensities were referred to the scale of the potentiometer. By suitable switching, 
the microammeter was used to check the current flowing in the potentiometer 
circuit. ‘This current, of 50 microamps., was found to change by less than 1 in 
500 during a night’s observations. As will be seen later, the method of observation 
served to eliminate the effect of any small uniform changes in this current except 
in so far as they might influence the extinction curve. Changes of I in 500 would 
be quite inappreciable here. The voltage supply for the multiplier (90 volts/stage) 
came from a miniature dry-battery pack which proved more suitable than a 
voltage-controlled rectifier supply used for earlier work. The batteries were 
satisfactory even after nearly two years of use. 
The linearity of response of the cell was tested in several fashions, viz. 
(i) by means of diaphragms of known transmission placed over the objective 
of a 6-inch refractor and using an actual star, 
(ii) by an inverse square method, 
(iii) by a summation method. 
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After a few preliminary trials method (i) was rejected partly because of atmospheric 
fluctuations and partly because of the difficulty of measuring the areas of the 
diaphragms to the degree of accuracy desired. Methods (ii) and (iii) necessitated 
the use of an artificial star. This consisted of an automobile lamp, mounted 
cap downwards and directed towards an inclined diffusing surface of white card 
from which light was scattered to a ground-glass screen and thence to a pin-hole. 
The ‘star’? was mounted on an improvised optical bench and the pin-hole was 
focused by means of a small lens on to a diffusing glass screen mounted in contact 
with the cell of the multiplier. For the method (ii) the “star”? was moved to 
definite distances (1 to 4 metres) from the focusing lens. The illumination of 
the lens and hence of the cell was then taken as being inversely proportional to 
the square of the distance. Various trials gave discrepancies of the order of 
o™-01 which was regarded as unsatisfactory. It was determined that these 
discrepancies arose in part from the alteration in size and position of the image 
of the star on the diffusing screen at the multiplier when the ‘star’? was moved 
into different positions. For the summation method (ii), the “star” was focused 
on the diffusing screen of the multiplier by a thin lens 10-5 cm. diameter and 80 cm. 
focallength. A blue filter was used here to avoid trouble from lack of achromatism 
of this lens. The lens could be covered by an opaque diaphragm in which had 
been drilled twelve holes each of diameter 14mm. These holes could be closed 
with light-proof stoppers when desired. ‘The test consisted in measuring the 
light admitted through the twelve holes simultaneously, then in two groups each 
of six different holes, in four groups each of three different holes, and finally 
through each of the twelve holes singly. If the cell was linear in response then 
the reading for all twelve holes simultaneously (W) was expected to equal the 
sum of the readings (G,, G, and G,) of each of the other arrangements. The 
following results were obtained as the means of several sets of readings when the 
illumination through twelve holes gave a multiplier current of 1-4 x 1077 amp. 
Two complete runs a and 6 were made for groups G, and G;. Thus any 
systematic deviation of linearity over a range of nearly three magnitudes would 
be less than 1 millimagnitude. 
TABLE | 
Group m W-—-G Mean 
(a) (b) (a and b) 
G, (6 hole) oM-75 o™ 9005 a o™ 9005 
G, (3 hole) 1-6 o™-9016 --0'0020 o™-9002 
Gs (1 hole) a9 —o™-9007 00012 o™ 9092 


(Column “ m”’ gives the ratio G/W expressed in magnitudes.) 


These results were not obtained without various precautions. Particular 
care was taken in the optical arrangements to ensure that the images through the 
various holes fell on the same part of the diffuser. ‘The lamp voltage was checked 
frequently, using a potentiometer reading to 0-0002 volts. ‘The lamp was run at 
5°8 volts, approximately one-third of this voltage being shunted to the potentio- 
meter for measuring. The accumulator supplying the lamp was charged 
continuously from a rectifier connected to the mains through a stabilizing trans- 
former so that the difference between the input and output currents did not 
exceed a few milliamps. This was necessary as the nature of the test demanded 
steady illumination over several hours. Further, to take account of any steady 
drift in the lamp output or the multiplier response, comparison observations were 
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suitably intercalated and a standard routine was adopted in exposing the cell and 
measuring the current. 

Throughout the programme the dark current was usually quite stable. 
After an exposure the meter needle would return to the zero position in the 
course of a minute or so and remain steady there. Occasionally, however, sudden 
bursts of dark current appeared and the cell then became unusable. This condi- 
tion was associated with warmer weather and the cell was restored to its usual 
behaviour by keeping it for a few hours in a domestic refrigerator. 

3. The programme.-—The programme aimed initially at measuring southern 
‘stars down to 2™-5 as well as a few northern stars. Actually some fainter southern 
stars were included. Known variable stars were omitted, as also were most 
double ‘stars, unless the components differed widely in magnitude. Double 
stars could give superposed spectra and render the wave-length limits meaningless. 
‘The star comparisons were divided into three groups termed primary, secondary 
and tertiary. 

The primary comparisons were measures made on seven stars forming a ring 
around the sky between — 20° and —70° declination. ‘They were restricted to 
stars showing the line H8 sufficiently well for setting. No double stars or stars 
reported as having variable radial velocity were included. Each primary star 
was compared with every other star in the primary group, so that 21 sets of 
comparisons were involved in this ring. ‘These seven stars were used as standards 
for the rest of the programme. 

The secondary group was composed of 24 stars in which the line Hf could 
usually be seen. Double stars for which the difference in brightness of the 
components exceeded 7™ were included, as errors due to the spectral uncertainty 
here would be small. Each of the secondary stars was measured against two 
primary stars, thus involving 48 sets of comparisons. 

The tertiary group included 32 stars in which Hf could not usually be seen 
and which therefore required setting by the fiducial mark. Also some double 
stars differing by less than 7™ in brightness were included here. Thirty-two sets 
of comparisons were used for this group. 

4. Methods of observation.—A single measurement of a star resulted from the 
mean of ten individual readings of the photometer made according to a regular 
plan over a period of 90 secs. ‘These individual readings showed some variation 
but this usually was less than +2 per cent on the class of night used. A complete 
measurement of a single star, including setting of the telescope, required 4 minutes 
with two observers: one guiding, the other reading the photometer. Each pair 
of stars was compared in the order a, 6, a to minimize the effect of any drift in the 
response of the apparatus or the transparéncy of the atmosphere. Comparisons 
were made at such times that the stars involved differed by less than 0-1 in sec s 
(z=zenith distance). This meant that corrections for atmospheric absorption 
generally were less than 0™-03 (see below). The order of each night’s observation 
was planned with the aid of a graph showing, for all programme stars, the variation 
of sec z with sidereal time. Each primary comparison was observed at least 10 
and usually 14 times spread over a minimum of 4 nights. Secondary and tertiary 
pairs were observed at least 6 times, and on at least 3 nights. Observations were 
made on 83 nights between 1948 September and 1950 April.* ‘To obtain 
extinction coefficients, nightly programmes were planned to include at least two 


* The telescope was available for this programme on half the nights between the dates mentioned. 
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stars, each of which could be measured at several different altitudes. The 
linearity and smoothness of the curves obtained by plotting log readings against 
sec z were criteria by which the night’s work was accepted or rejected. 

5. Filter and diaphragm calibrations.—Although in the later part of the 
programme the observations were arranged so that a knowledge of the optical 
densities of the filters and diaphragm was not needed, it was necessary to determine 
these quantities for some of the earlier work. Densities were involved in only 
18 comparisons (out of a total of 101). Even in these comparisons 39 per ceni 
of the observations were reduced without the need of optical densities. 

Measurements of the filter density showed that from 1948 September to 
late in 1949 April the o™-5 filter maintained a constant density (0™-587 + 0-001 s.d., 
10 observations). Between 1949 April 23 and 1949 April 26 it decreased suddenly 
in density and then remained effectively constant again (at 0™-556 + 0-004 s.d., 
8 observations) till the completion of the observations. The determinations of 
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Fic. 3.— Absorption Curves for a Typical Night. 


the optical density of the 2™ filter exhibited rather irregular scatter and a mean 
value (of 1™-o79 + 0-009 s.d., 12 observations) was adopted throughout. As has 
been described above, the transmission factor of the Maltese-cross diaphragm 
could be expected to change for different series of observations. Accordingly 
individual values were adopted for the density of this diaphragm for the various 
series. These values ranged from 0™-898 to 0™-g26. Only eight comparisons 
involved the Maltese-cross diaphragm. 

6. Atmospheric extinction.—Absorption curves were obtained on 67 nights 
out of the 76 used in the final reduction. On most nights two or more curves 
were available. Curves for a typical night are shown in Fig. 3. ‘The mean 
extinction as determined with rising stars differed slightly from that with setting 
stars. The means for 26 nights (i.e. all nights for which eastern and western 
curves were both available) gave the apparent rising extinction as 0-254 and setting 


23 
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extinction as 0-281 (unit = magnitude/sec z). The difference, 0-027 + 0-009 (s.d.) 
was statistically significant. 
Several reasons might be advanced for the existence of this difference, viz. 
(i) a true difference in absorption in eastern and western skies, 
(ii) a fatigue effect in the cell response or some other instrumental drift, 
(iii) an apparent difference due to one or both of the absorptions increasing 
with time, 


(iv) an absorption varying with zenith distance differently in east and west. 


There are no obviously important sources of local pollution to account for (i) 
nor is it thought that fatigue effects are significant with the low light-levels and 
short exposures used. In examining the results it was not practicable to separate 
these effects. Four methods of examination however presented themselves, viz. 


(A) To adopt for each night the mean of all the coefficients measured on that 
night. 


(B) To adopt for all nights the grand mean of all the coefficients measured 
in the programme. 


(C) To adopt separate means derived from rising and setting stars (e. g. 0°254 
and 0-281). 


(D) To adopt coefficients adjusted to take account of variations of absorption 
coefficient with secz. For this purpose it was assumed that the smaller 
mean extinction coefficient measured with rising stars (k,) represented 
a value which was constant for all stars before transit. For stars after 
transit it was assumed that the absorption coefficient for a particular 
secz was k,=k,+6(secz—1). Values of k, and 5 derived from the 
observations were 0-254 and 0-014 respectively. This value of 56 was 
insufficient to bring about any detectable departure from linearity in 
the absorption curve. 

Some of these methods, because of the distribution of the observations, take 
account of more than one of the effects listed as (i)-(iv). Thus Method A (if rising 
and setting stars were equally represented in the mean) would eliminate effects (ii) 
and (iii) from the observations as a whole. Method B would similarly take account 
of effects (ii) and (iii) and might be expected to give smoother results if the 
determinations of the absorption coefficients were greatly affected by individual 
errors. Method C might take acceunt of effect (i) but was rejected because it 
involved an abrupt change at the zenith. Method D would approximately 
account for one form of effects (iii) and (iv). Methods A and D were applied 
to the 21 comparisons of primary stars. It was found that the average deviation 
in the mean magnitudes derived by the two methods was 0™-005, no single 
deviation exceeding 0™-o10. The residuals of observed comparisons relative to 
values expected from assigned magnitudes however were rather less for Method A 
than Method D. Further, Method A gave results showing less scatter than did 
Method B. Accordingly Method A was used to reduce all the results except 
for nine nights when no absorption coefficient was determined and the grand 
mean coefficient (B) was used instead. 

7. Reductions and results.—After preliminary reduction certain measurements 
were rejected either because of cloud or haziness noted at the time of measurement 
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or because of irregularities in the extinction curve. Out of about 2500 measure- 
ments some 250 were not used for these reasons. Most of the rejected measures 
were made on seven unsatisfactory nights. All the readings taken on these 
nights were rejected. 

Magnitudes for the seven primary stars were obtained from a least squares 
solution for the 21 comparisons, an arbitrary zero being adopted for the magnitude 
scale. The solution aimed at assigning a set of magnitudes to the seven stars 
which would best represent the 21 observed magnitude differences. Next a 
magnitude was obtained for each secondary star as the mean of the results of the 
comparisons with the two primary stars. Each tertiary magnitude was found 
from the single comparison with a primary star. The internal uncertainty 
attaching to these results was estimated as follows : 


(i) Primary stars s.d.= +0™-002. 
(ii) Secondary stars = +0™-005,. 
(iii) ‘Tertiary stars = +0™-008. 


It was desired to link the present observations with the zero of the North 
Polar sequence. ‘Two methods presented themselves, viz. 


(a) by comparison with the Harvard Mimeogram magnitudes, 
(6) by comparison with certain results given by Eggen. 


The Harvard Mimeogram results (4) were brought to the International Photo- 
visual Scale (IPv) by the equation given by C. Payne-Gaposchkin, 


IPv(h) = HPv + 0-10 + 0-06C,. (1) 


The colour index used in this equation (C,=1Pv—IPg) was here computed for 
each star from the spectrophotometric gradients G given by Gascoigne (2) or 
else by Greenwich (§). The definition of G gives 


IPv —IPg =1-086(1/A,,, — 1/A,,)G. 


The effective wave-lengths A,, and A,,. were taken from the work of Seares and 
Joyner (6). It was known that colour index was not strictly proportional to 
gradient, but using the above computed values of C; in equation (1) did not 
introduce any appreciable error as the colour correction was so small. 

The photoelectric results (m4;59) were also reduced to the wave-length A,,, by 
using the expression 


IPv(s) = my559 + 0°25 + 1°086(1/A,,. — 1/4550)G. 


Although the gradients employed were determined with band widths of only 
2-3 A. they had been shown by Gascoigne to agree well with the gradients measured 
photoelectrically by J. S. Hall (7) using a band width of 4004. Presumably no 
error would be made in applying them to the 100 a. band used here. 

Fig. 4(a) compares the I Pv values derived from both the Mimeogram and the 
M4559 results. The differences IPv(h)—IPv(s) were plotted against IPv(s) for 
individual stars. Double stars with components differing in brightness by less 
than 5™ have been marked with a cross. The values IPv(h)—IPv(s) were also 
plotted against gradient (Fig. 5). The results showed that whilst there was no 
systematic colour effect, the two magnitude scales deviated at the bright end. 
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Fic. 4.—Comparison with Harvard, Stromlo Photographic and Lick Results. 
(a) IPv(h)—IPv(s) vs. IP v(s). 


(5) Mgsso(PB)— Mgsso(PE) VS. ™y;5(Pe). 
(c) Myggo( Lick) —m,,,9(Stromlo) vs. m,,,,(Stromlo). 


Least squares solutions gave the constants in the following equations :— 
(i) IPv(A) = IPv(s) — 0-09 + 9-111 Pv(s) — 0-025(1Pv(s))*, 


0-126 + 0°029 (s.d) 


(ii) IPv(h) = 1Pv(s) + 0-04 — antilog 0-41Pv(s) * 
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The second equation gave a slightly better fit than the first, the s.d. of one residual 
being (i) +0™-o151 and (ii) +0™-0146. The second equation was derived by 
assuming that the actual light intensities after correction for colour (L, and L,) 
were related by the expression 

L,=aLp. 
This is of course a purely empirical relation. 

As the linearity of the photometer scale had been carefully examined it was 
unlikely that material scale error existed in the photoelectric results. This view 
was upheld by a comparison of the m4,» results with some unpublished measures 
generously made available by Professor Woolley. ‘These results were obtained 
by an improved version of the method described by Woolley and Gottlieb (8). 
They have been determined photographically with a band width of 250a. Their 
scale was established with the aid of carefully measured sectors, and colours were 
obtained. The two scales were compared at A=4550a. and the difference 
(photographic—photoelectric) was plotted against the photoelectric results (see 
Fig. 46). There was no systematic deviation with magnitude and the mean 
deviation over a range of 4™ was + 0™-04. 
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Fic. 5.—IPv(h)—IPv(s) vs. Gradient. 


The m4;59 scale also agreed with the above-mentioned photoelectric scale 
given by Eggen. Eggen’s results were expressed on a system (Pg,) which was 
directly linked with the North Polar Sequence (9) so that the zero-points agreed 
between the sixth and eighth magnitudes. 

The results for the eight stars available for comparison (10) gave the following 
relation :-— 


M4559 = Pg, —0°16C,,. 
In this equation C, is a colour defined by Eggen. 

Table II shows values of Pg, and C, as well as figures for m4559 as computed 
from Eggen’s results and as observed. The deviations between the two sets of 
results (computed—observed) are plotted in Fig. 4(c). ‘The maximum deviation 
is O™-03 over a range of 5™ in brightness and 1™ in colour. 
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Because of this agreement the zero-point for the m,;59 catalogue (‘l'able III) 
has been chosen to agree with the zero-point of the Pg, scale, i.e. the M4559 results, 
when corrected for colour would be expected to agree with the I Pg scale at the 
mean magnitude of the six stars of the North Polar Sequence (7™) chosen by 
Eggen to define the zero of his scale. ‘This means that the listed m,;59 figures 
are in general 0™-25 brighter than magnitudes which could be derived by reducing 
the Harvard Pv system to A,~=4550. The most notable exception to this 
statement is the magnitude of « C Maj which on the photoelectric scale is 
™4559 = — 1°49 whilst the derived figure is — 1°53. Deviations of over o™-1 have 
been noted between photographic and photoelectric results for the brighter 
stars of the North Polar Sequence (9, 11, 12). 


TABLE II 
Values of mys59 as observed and as computed from Eggen’s results 
Star Obs. Pgp Cp M4550 
computed observed 
y Gemi 1°88 —0'065 1°89 1°86 
a C Maj s —1°53 —0'07 —1°52 — 1°49 
B Leon ; 2°11 +o0°08 "10 2°10 
a Cor B s 2°12 —0°05 "33 ‘II 
a Ophi 2°11 +o'ro 09 ‘10 
a Aqil 0°88 +0°165 0°85 83 
€ Pegs 3°74 + 1°37 3°52 "53 
a Psc A p I'll +O°O15 ess ‘II 


8. Catalogue—Table III gives the magnitudes at A=4550 for 63 stars. 
The magnitudes determined by primary, secondary and tertiary methods are 


respectively denoted by the letters p, s and t. Double stars, the components 
of which differ in brightness by less than 5™, are denoted by (d) as it is considered 
that the spectral range for these results may be outside the selected band. 


Taste III 

Catalogue 

Star a (1950) 5 (1950) 

hm <a 

a Andr 00 06 +28 49 

B Hydi 00 23 77 32 

a Phoe 00 24 42 35 
B Ceti 00 41 —18 16 708 
« Erid o1 36 57 29 ‘274 
x Hydi OI 57 61 49 ‘998 
y Hydi 03 48 74 24 *369 
B Orio 05 12 o8 15 138 
y Orio o5 22 +o6 18 434. 
¢ Orio o5 38 or 58 598 
« Orio ©5 45 09 41 *g22 
B C Maj 20 17 56 “792 
a Cari 23 52 40 —o'686 
y Gemi 06 35 27 857 
a C Maj 06 43 39 488 
a Pict 06 48 53 3°326 
e C Maj 06 57 3 54 295 
8 C Maj 07 06 19 "305 
» C Maj 07 22 12 347 
a C Min ©7 37 5 21 608 
{ Pupp 08 02 - 52 032 
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TABLE III (cont.) 
Catalogue 
Star x (1950) ™M4550 


m 
y Velr 08 —47 1°571 
€ Cari 21 —59 2°794 
5 Velr 43 —54 “B80 
A Velr 06 —43 "395 
fb Cari 13 —69 596 
t Cari 16 —59 318 
« Velr 21 -54 ‘290 
x Hyda 25 08 ‘002 
x Leon 06 12 238 
~ Leon 47 -14 § 101 
a Cruc 24 —62 536 
y Cruc 28 56 s 642 
y Cent 39 —48 038 
6 Cruc 45 -59 035 
a Virg 4-23 —10 § ‘773 
p Cent oo 60 "402 
6 Cent 04 -36 ‘716 
« Boot 13 -19 ‘780 
a Cent 36 -~60 3 114 
y Tr Au 5 14 ~68 ‘800 
a Cor B 33 26 53 114 
6 Scor < <7 22 2 ‘202 
a Scor 26 26 323 
x Tr Au 43 68 : 903 
€ Scor 47 — 34 032 
6 Arae 27 60 3 454 
A Scor 30 —37 *380 
a Ophi 33 12 3 -096 
6 Scor 34 —42 5 ‘081 
« Scor 39 ~39 +196 
€ Sgtr 21 -34 2 ‘742 
o Sgtr 52 —26 883 
« Aqil 48 -o8 834 
5 Pavo —66 ‘O75 
a Pavo 56 5 “718 
f Pavo 66 "451 
e Pegs 21 Log + "531 
x Grus 22 O8 47 585 
6 Grus -47 "363 
B Octn 22 41 —81 3 ‘179 
a Psc A 22 55 29 53 ‘108 
« Pegs 23 +14 56 "405 

* Magnitudes based on three observations only. 

+ These stars are listed by Giissow and Guthnick (13) as having slight variation of 
amplitude of the order of o™-1. In the present work, the standard deviations of 
individual comparisons for aC Min and a Leon were amongst the highest found 
(o™-036 and o™-o41 respectively). The stars « Cor B and «a Virg showed only normal 
scatter (002 s.d. for individual comparisons). As, in general, such comparisons were 
made on from three to six nights, variations could escape detection. 
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SPECTROSCOPIC BINARIES 
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Summary 


Observations of six southern spectroscopic binary stars were made with 
an electron multiplier photometer. Of the stars examined a Cari, 6* Cruc 
and vy Cent were found to have no regular light variation with the spectro- 
scopic period, 6* Cruc shows some irregular variation of amplitude +0™-o2. 
66 Erid and 19 Erid were not fully covered but no variations were detected. 
ao Scor was found to have a variation similar to that of a dwarf Cepheid. 
yp Cent (which is not actually a spectroscopic binary but was used as a 
comparison star) was found to vary in an irregular fashion with a range of 
om-s. 





Suitable southern spectroscopic binaries were selected from Moore’s Fourth 
Catalogue (1) for photoelectric observation in the expectation of discovering 
shallow eclipses. ‘The programme included the stars set out in Table I which also 
gives the comparison stars used. ‘The comparison stars were selected to be, as 
far as possible, of the same spectral class as the variable so as to avoid differential 
absorption effects. The observations were made almost entirely with an electron 


multiplier photometer attached to the Farnham 6-inch refractor. This instrument 
and the technique adopted have been described in a previous paper dealing with 
the eclipsing variable { Phoe (2) which had also been included in the present 
programme. A few observations were made with a multiplier photometer on 
the Oddie 9-inch refractor (3). 


TABLE | 
Binary and comparison stars 
(1900 Positions) 


Star R.A. Dec. Mag. Sp. Period Comparison R.A. Dec. Mag. 5 

19 Erid “ Bi 21 58 4°32 B8 "227 17 Erid of 26 5 25 4°80 

66 Erid 05 o2 447 5'19 Bo 5°5 62 Erid 04 51 5 5°46 

a Cari* 09 08 58 33 3°56 B3 . 36631 og 09 661 4°18 

62? Cruc 11 59 62 36 4°98 B3 3° 4573T I! 53 61 53 5°70 

v Cent 12 36 ar 1% 3°33; Ba 2°62 ¢ Cent ry 33 41 37 4°05 

o Scor 16 15 25 21 3°08 Br 2 7 Scor 16 30 —28 2°91 

p Cent 13 44 —41 59 3°32 Bap... ¢ Cent 13 52 —42 15 3°53 
Fig. I gives results for some of these stars. No systematic variation is shown 

in 66 Erid, cf. J. Stebbins (4), 19 Erid, a Cari, #* Cruc or v Cent, cf. R. E. Wilson 

(5). Comparison of the curves for #* Cruc and v Cent, which were obtained by 

the same technique, suggests that as the former shows appreciably more scatter, 

there may be some irregular small-scale variations in the light of either 6* Cruc 

or its comparison star. No definite periodicity could be determined. 


* L. Boss, P.G.C., 2473, 1910. 
+ F. Schlesinger, Catalogue of Bright Stars, 1930. 
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o Scor.—This system exhibits two types of radial velocity variation, viz. a 
change with the very short period of 0°-2468 and a change with a longer and some- 
what varying period of about 33° (Henroteau (6)). Strictly speaking the short- 
period variation dealt with here is most likely not of the binary type and the star 
is included in Moore’s Catalogue because of the long-period variation, which 
however has not been examined photometrically. Henroteau classifies the star 
as being of the 8 C Maj type on account of its spectroscopic changes. Such stars 
are characterized by a period of a few hours and light variations of very small 
amplitude have been detected in some members of the class. Stebbins (4) has 
compared the variations of o Scor with 7 Scor and concluded that there was some 
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Fic. 1.—The light curves of some spectroscopic binaries. The ordinates represent the 
magnitude differences between the binary and the comparison star listed in Table I. The figures 
on the right show the average number of comparisons meaned to obtain each of the plotted points. 
The curve marked a Scor (B-Y) gives the colour change of the star. 


kind of variation with a range of at least o™-03 but the cbservations were insufficient 
to derive a complete light curve. Three sets of measures have been made in 
the present case, viz. :— 
(a) A series on seven nights in 1947 using a 931A multiplier on the 6-inch 
Farnham refractor without any filter on the cell. 
(b) A series on two nights in 1950 using a 1P21 multiplier on the g-inch 
Oddie refractor in two colours corresponding to a blue filter (Schott 
F3653, 2mm.) and a yellow filter (Schott F4937, 2 mm.). 
The results of all sets showed a variation which could be represented by a period 
of o*-246841 derived from the epochs of minimum light for the 1947 and 1950 
observations, viz. :— 


J.D. =2432 28408 and J.D.=2433 396-10, 
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by assuming that 4505 cycles elapsed in that interval. ‘The spectroscopic period 
given by Henroteau, viz. 0'-246834 required 4505-2 cycles. Owing to the 
considerable lapse of time between the spectroscopic and photometric observations 
and the very short period it is inadvisable to attempt to identify the two sets of 
results. 

The earlier series is shown in Fig. 1, by the curve marked “cell”. The 
normals are the means of about twenty comparisons, with standard deviations 
about 0-006 per normal. The range is o™-08 and the form of the light variation 
is similar to that shown by other 8 C Maj stars. The two-colour curves are 
also shown in Fig. 1 marked “yellow” and “blue” respectively; here the points 
are means of four comparisons. ‘These confirm the earlier results. The curve 
marked B-Y shows the colour. Corrections were applied to bring the colour 
to the zenith. The curve makes the star bluer at maximum light as is the case 
for Cepheids generally. 

» Cent.—This star exhibits variable line emission (for references see Merrill 
and Burwell (7)). Early in the present century it showed the Balmer series in 
-emission as far as He and also bright lines of ionized iron, but by 1918 all bright 
hydrogen had disappeared. Gerasimovic (8) investigated the photographic 
density of plates exposed at times when the bright lines were present and absent. 
He concluded from measures on six plates that the star would be fainter by 
‘oI or o"-2 when the lines were not visible. 
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Fic. 2.—Light variations of 4 Cent. The points for the upper curves are the means of about 
25 comparisons. The lowest diagram marked p Cent—p Cent represents the mutual comparison of 
successive observations of ~ Cent. Thus measures a, b and : of Cent have been compared in 
the ratio b: 4a+c). The deviations of this comparison from zero magnitude show the instru- 
mental and atmospheric scatter which, for the night plotted, is represented by a standard deviation 
of o™-006 per comparison. This value is below average. 








In the Spring of 1947 irregularities were noticed in the comparison of v Cent 
and » Cent. In the light of later information these irregularities were attributed 
to changes in the light from yz Cent. No short period could be discovered, although 
irregular variations amounting to 0-05 were found to occur in the course of a few 
days or less. Observations were carried out for the 1948 and 1949 seasons with 
¢ Cent as a comparison star and again no short-periodic variations of a regular 
nature could be found. The magnitude difference between v Cent and ¢ Cent 
was found to be constant (0-428) and this figure was used to bring all results 
toacommon basis. ‘The observations are shown plotted against the time in Fig. 2, 
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where the normals are means of about 25 comparisons. ‘There is a slow variation 
which in the 1948 results seems to have a period of about 60 days but this is not 
persistent and the general picture is of an irregular variable. The light level for 
the 1949 season was about 0"-1 above that for the previous year and the complete 
range of variation was 0-2. Unfortunately the complete period of the obser- 
vations was not covered by spectrograms. Dr Gascoigne has kindly made 
available some objective prism spectrograms taken on J.D. 2,432,754; 2,433,012; 
2,433,020. These show a bright line superposed on the absorption line 
at HB. ‘The emission line is displaced somewhat to the red of the centre of the 
broad absorption band. There is no difference in the three spectrograms 
available that might not be attributed to seeing conditions and it is not possible 
to make any statement about a correlation between the light output and the 
spectral changes except that the o"-1 change of light level between 1948 and 1949 
was not associated with any marked changes in the emission lines. 


Commonwealth Observatory, 
Mount Stromlo, 
Canberra, Australia : 
1950 December 20. 
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OCCULTATIONS OF STARS BY THE MOON OBSERVED AT 
THE NIZAMIAH OBSERVATORY, HYDERABAD, 
DURING THE YEAR 1950 


(Communicated by the Director) 
(Received 1951 February 9) 


TABLE I 
Phase Date U.T. 
1950 
h m 8 
Jan. 33 22 20°4 
2 13 2°7 
29 14 21°9 
29 18 5'2 
29 18 38°7 
31 21 55°! 
4 16 4°38 
5 23 25 36% 
6 ° 49°00 
2 13 59°5 
24 13 4471 
24 13 56°2 
24 14 55°6 
24 14 °0°6 
2 14 2°5 
28 14 6°9 
I 14 39°6 
8 23 43 52°9 
I 18 57°0 
2 14 31°5 
28 14 52°! 
22 14 34°0 
19 14 § 6-7 
13 18°8 
20 17 37°8 M.G. 
20 18 o'8 M.G. 
24 13 53°5 M.G. 
15 57°0 M.G. 
16 15 14°8 M.G. 
18 14 55°2 M.G. 
20 12 32°5 M.G. 
20 17 14 11°4 M.G. 
20 17 47 50°5 M.G. 
21 14 26 20°'8 M.G. 
*% 21 14 34 §2°0 M.G. 
és 21 14 38 25°8 M.G. 


560 

562 

° §62 

| §67 

* 564 

578 

1149 

1279 

2157 

1712 

1181 

1644 

2538 

2677 

3225 
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3240 

146 

3405 

3520 

203 

425 

440 

N.Z.C. 438 
N.Z.C. 541 
N.Z.C. 542 
N.Z.C. 543 
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The last column indicates the initials of the observer, Mahomed Ghouse. 


The instrument used throughout was the 15-inch Grubb refractor. The 
corrections to the standard clock were determined from observations with the 
transit instrument. The reductions were carried out by S. Aravamuthan 
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assisted by V. V. Raghaviah on a Facit calculating machine, after suitably modi- 
fying the formulae given in the Supplement to the Nautical Almanac for 1938. 
A correction of —0"-00076 was applied to the observed time, to allow for the error 
of —1”-5 in the mean longitude of the Moon. 


Tasze II 
cos (p—) sin (p—) No. cos (p—y) sin (p—) 


Z 


—o'66 +-0°75 19 1° 1°00 —0'03 
+ O°47 +088 20 0°76 + 0°65 
+-0°99 +o°12 21 1°00 + 0°03 
+0°58 -o81 22 +O°34 ~0'94 

0°16 —0'99 23 0°95 —0°30 
+096 —0'28 24 + 0°96 +0°27 

o'92 0°39 25 +080 -0°59 
}-0°97 0°24 26 0°98 o'18 

0°76 + 0°65 27 + 1°00 -0'O3 
+-0°99 —O°Is 28 0°96 +-0°27 
+ 0°97 —0'23 29 +-O°43 +o'g0 
+ 0°62 0°78 30 +-0'87 0°49 
+-o'84 O'54 31 1°00 0°04 
+ 0°44 o'go 32 0°99 —o'16 
+ 0°99 —Or'll 33 +-O°40 + 0°92 
+ 0°97 — 0°25 34 0°89 ~O°45 

0°49 0°87 35 0°98 +or21 

0°79 —o'61 36 1°00 0°05 
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The mean places for 1950-0 of the five stars (Nos. 10, 12, 13, 14 and 15) are 


bh J 
ae oe 
= + 


given in Table III. For star Nos. 10, 12, 14 and 15, proper motions have been 
applied to the coordinates given in the catalogue in computing the mean place 
for use in the reduction. 


Taste III 
R.A. (1950°0) Dec. (1950°0) Remarks 
h m ~ es » ” 
3 46 17°68 +24 14 42°5 A. G. Berlin B1201, 
Hy +08-0017 ; Ug—0” 040 
(Yale) 
3 46 4 05 48°7 A. G. Berlin B1203, 
Hg + 080015 : fg—0"'038 
(Yale) 
3 46 Boss GC4609 
A. G. Berlin B1208, 
Na+ 08-0003 ; [Wg+0" 010 
(Yale) 
A. G. Berlin B1218, 
My +o8-0015 ; @g—O” O12 
(Yale) 


Nizamiah Observatory, 
Hyderabad, 
Deccan, India: 
1951 February 1. 
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